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Abstract: The isospin-breaking decay r] — > 3tt is an ideal tool to extract information on 
light quark mass ratios from experiment. For a precise determination, however, a detailed 
description of the Dalitz plot distribution is necessary. In that respect, in particular the 
slope parameter a of the neutral decay channel causes some concern, since the one-loop 
prediction from chiral perturbation theory misses the experimental value substantially. We 
use the modified non-relativistic effective field-theory, a dedicated framework to analyze 
final-state interactions beyond one loop including isospin-breaking corrections, to extract 
charged and neutral Dalitz plot parameters. Matching to chiral perturbation theory at 
next-to-leading order, we find a = —0.025 ± 0.005, in marginal agreement with experi- 
mental findings. We derive a relation between charged and neutral decay parameters that 
points towards a significant tension between the most recent KLOE measurements of these 
observables. 
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1. Introduction 



The decay rj — > 3n has been the center of attention in many theoretical and experimental 
works over the recent decades. The considerable interest is due to the fact that the decay 
can only occur via isospin-breaking operators and is therefore sensitive to the up- and down- 
quark mass difference. Indeed, the rj — > 3tt transition amplitude is inversely proportional 
to the quark mass double ratio Q 2 , 

1 m 2 , - m\ 1 

7t2 = ~~ 2 ^2 ' rn = -{m u + m d ), (1.1) 

£r raj. — m z 2 

and thus the decay provides an excellent testing ground for the breaking of chiral symmetry. 

Despite valiant efforts it seemed difficult to bring theoretical description and exper- 
imental results in agreement. First attempts that relied on an electromagnetic transi- 
tion [1,2] were unsuccessful in explaining the decay. SU(3) current algebra techniques in 
combination with the partially conserved axial- vector current hypothesis [3,4] were gener- 
alized to SU(3) chiral perturbation theory (ChPT) and initiated systematic improvements 
to the decay rate. While the one- and two-loop corrections to the decay were sizable [5,6], a 
consistent implementation of electromagnetic contributions only lead to small effects [7,8]. 
Despite these theoretical improvements the Dalitz plot expansion, especially of the neutral 
decay, remained an unsolved puzzle. The slope a vanishes at leading order, while at next- 
to-leading order (0(p 4 ), one loop) it disagrees in sign with experimental findings [9-17]. 
The same holds for the next-to-next-to-leading order (0(p 6 ), two loops) calculation [6]. 
The error on the final result is rather large, so that it allows for a negative slope parame- 
ter. However, this error is not based on the uncertainties due to the low-energy constants 
at 0(p 6 ), which are estimated by resonance saturation, but results solely from the authors' 
fitting procedure. 

It has been argued that tttt final-state interactions are the dominant force behind the 
sizable corrections [18,19], motivating several dispersive analyses [20-23] (see also Ref. [24]), 
which were able to give a more robust prediction of the slope parameter. Among the short- 
comings of these dispersion relation techniques and the next-to-next-to-leading-order cal- 
culation is the treatment of higher-order isospin-breaking effects due to electromagnetism, 
as for example the mass difference between charged and neutral pions. It is not yet clear 
how to incorporate these effects. 

An analysis of rj — > 3-ir in the framework of unitarized chiral perturbation theory 
has been conducted in Ref. [25], producing remarkable agreement with experiment. In 
particular, the experimental value of the slope parameter in the neutral decay channel can 
be accommodated. However, since this approach is based on an elaborate fitting procedure, 
wherein the U(3) expansion parameters are determined from several hadronic rj and rj' decay 
channels, among those rj — > Sir, we do not consider this value for the slope parameter an 
unbiased prediction. Finally, a study of rj — > Sir in the framework of resummed ChPT is 
currently work in progress [26]. 

In this work we attempt to bridge the gap between the ChPT prediction and the disper- 
sive analysis using the modified non-relativistic effective field theory framework (NREFT) . 
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While this framework does not allow for a fundamental prediction of physical observables, it 
is ideally suited to study the dynamics of the final-state interactions. At two-loop accuracy 
and with the correct empirical tttt scattering parameters, we ought to have a reasonable 
approximation to the full dispersive resummation of rescattering effects at hand, so when 
matching to ChPT at 0{p A ), we can hope to find a transparent interpretation of the disper- 
sive results obtained in a similar fashion [20]. Additionally, the non-relativistic framework 
provides access to investigating the effects of isospin-breaking corrections. 

This article is organized as follows. In Sect. 2 we begin with a short description of the 
Dalitz plot expansion and the conventions used throughout this work. An introduction to 
the non-relativistic framework, its power counting, the matching procedure, and numerical 
input is given in Sect. 3. Sections 4 and 5 comprise the analytic and numerical results in the 
isospin limit and with isospin breaking included. In Sect. 6 we study final-state interaction 
effects on an isospin relation between the charged and the neutral decay channel. In 
Sect. 7 finally, we briefly comment on the 77 — > 37r partial widths and their ratio, before 
summarizing our findings in Sect. 8. Several of the more laborious formulae are relegated 
to the appendices. 



2. Dalitz plot expansion of the decay amplitude 

In the following we consider the charged and neutral decay modes 

rj(P v ) -»• Tr + (pi)ir~(p 2 )ir°(p 3 ) , Si + s 2 + s 3 = 3s c = M 2 + 2M* + M* , 

Q C = M V - 2M W - M^o , 

r){P v ) -> Tr ( Pl )ir (p 2 )iT (p 3 ) , Sl + s 2 + s 3 = 3s„ = M 2 + 3M 2 , 

Q n = M v - 3M n0 , (2.1) 

where the kinematical variables are defined as s, = {P-q—Pi) 2 with pf = Mf, i = 1,2, 3, and 
Q n /c is the excess energy of the respective channel. 1 We will use the notation M w = M 7r ± 
throughout. 

In experimental analyses of these decays, the squared absolute value of the amplitude 
is conventionally expanded as a polynomial around the center of the Dalitz plot in terms 
of symmetrized coordinates. For the charged decay channel one uses 

r-E 2 -Ei si - s 2 3E 3 s n -s 3 

X = ^^"73X' y = ^;- 1 = — +6 ' (2 - 2) 

where Ei = p® — Mi is the kinetic energy of i-th particle in the rj rest frame, and we used 
the definitions 

n Ml + Mf - Si 2 Q n . 1 

= _ji 1 1 R = _ Mri Q c/n , 5=^-1 . 2.3 



1 For convenience, we use a different notation from what is usually found in the literature. The transition 
can be made setting s\ — t, S2 = u, S3 = s. 
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For the neutral channel one defines 



2,2 
x n ~+~ Vn i 



Vicos(0) , y n = i/zsin((f)) , 



(2.4) 



where we have introduced polar coordinates in the center of the Dalitz plot. These defi- 
nitions of x n and y n agree with x and y only for M v = M n o. Experimental data is then 
fitted to the Dalitz plot distribution, which is of the form (assuming charge conjugation 
invariance) 



\M c {x,y)\ 2 = \M c \ 2 {l+ay + by 2 + dx 2 + fy 3 + gx 2 y + . . .} , 
\M n (z)\ 2 = \AT n \ 2 {l + 2az + 2/3z 3 / 2 sin(30) + 2jz 2 + ...}, 



(2.5) 



1.0 r 



where a, b, d, /, g and a, /3, 7 are the Dalitz plot parameters and jV c , M n are the normaliza- 
tions of the charged and the neutral decay, respectively. We note that of the higher-order 
parameters beyond quadratic order in x and y, only / has been measured so far (by the 
KLOE collaboration [27]). However, with the advent of very high statistics measurements 
for rj —> 37T° e.g. at MAMI [28], a determination of (3 and 7 might not be beyond the realm 
of possibility. 

We wish to comment on the validity 
of the polynomial expansion Eq. (2.5) 
in particular for the neutral decay chan- 
nel. The boundary of the Dalitz plot for 
rj — > 3tt° is shown in Fig. 1. The dotted 
lines denote the three symmetry axes, 
the dotted circle depicts the beginning 
of the rapid decrease of pure phase space 
for radii y/z > \/0.756. It is important 
to note that the cusps due to ir + ir~ —> 
7r°7r° final state rescattering occur at 
Si = 4M 2 and not at a single z value; 
the smallest and the largest values of z 
crossing the cusp lines (z = 0.598 and 
z = 0.882, respectively) are indicated 
at the corresponding arrows. Therefore 
the polynomial representation for the 
neutral Dalitz plot distribution (2.5) is 
only valid for z < 0.598, i.e. inside the 
dashed circle. 

Table 1 summarizes the latest experimental determinations and theoretical predictions 
for a. In the following we propose an explanation for the disagreement between the ChPT 
result and experimental data. Our findings substantiate the dispersive result [20], and 




Figure 1: Boundary of the r\ — > 3ir° Dalitz plot. 
Dotted: symmetry axes and biggest enclosed circle. 
Dashed: cusps at Sj = 4M^ and corresponding cir- 
cle. Arrows: indicating specific z values (see text for 
details). 
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we are confident that it leads to a better understanding of the nature of the final-state 
interactions. 

It is worthwhile at this point to quote the ChPT decay amplitudes at leading order p 2 
and up to next-to- leading order in the isospin- breaking parameters — m u and e 2 . For 
the charged and neutral decay, respectively, they read [8] (we use the Condon-Shortley 
phase convention throughout) 



where F n = 92.2 MeV is the pion decay constant, and Bq is linked to the quark condensate 
in the (SU(3)) chiral limit in the standard manner. Equation (2.6) shows the isospin- 
violating nature of the decay, as both leading-order amplitudes are explicitly of order 
m-d — m u- At that order in isospin breaking, the rj — > 3-ir amplitudes fulfill the well-known 
AI = 1 relation 



which can be easily checked in Eq. (2.6). This relation even holds in general at leading 
order in the isospin-breaking parameters, i.e. also for terms of 0(e 2 ) [7], and is only violated 
at 0((md — m u )e 2 ) [8]. In the following, we will often adopt a loose way of talking and 
speak about the isospin limit for the charged and neutral r\ — > Sir amplitudes; this only 
refers to the approximation in which the relation Eq. (2.7) holds, in particular M n o = M n , 
and not to the limit m u = where the decay n — > Sir is (almost) forbidden. 

Note furthermore that all contributions involving = M 2 — M 2 = 0(e 2 ) in the 
charged decay amplitude have been absorbed by writing Eq. (2.6) in terms of s n . This 
motivates an expansion of the decay amplitudes of both channels around the point S3 = s n , 
si = S2'- we anticipate that, defined this way, higher-order isospin-breaking corrections to 
the AI = 1 rule for the normalization of the amplitude are going to be of chiral order p , 
without contributions from the tree-level amplitudes Eq. (2.6), and therefore small. This 
"center" of the Dalitz plot then corresponds to s± = S2 = S3 = s n and x n = y n = z = 
in the neutral channel, but to s\ = S2 = s n + A n , s 3 = s n or x = and y = 5 / in the 
charged case. The charged and neutral decay amplitudes then take the form 




M^°(S 1 ,S2,S 3 ) 




(2.6) 



M n (si,s 2 , s 3 ) = -M c (s 1 ,s 2 , s 3 ) - M c (s 2 , s 3 , si) - -M c (s 3 , si, s 2 ) 



(2.7) 




S2? + f(s 3 - Sn) 3 




+ 0((R n -R c ) 2 )} , 

Snf] 
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Theory 




a 




ChPT 0{p 4 ) [5] 


+0.013 






ChPT 0(p 6 ) [6] 


+0.013 


± 


0.032 


Dispersive [20] 


-0.007 




-0.014 


Experiment 




a 




Crystal Ball@BNL [9] 


-0.031 


± 


0.004 


Crystal Barrel@LEAR [10] 


-0.052 


± 


0.020 


GAMS-2000 [11] 


-0.022 


± 


0.023 


KLOE [12] 


-0.0301 


± 


0-0035l° : 0022 


A/T A 1\/TT T3 ri Ql 

MA1V11-D [Lo\ 


U.UoZ 


1 

± 


U.UUz ± U.UUz 


MAMI-C [14] 


0.032 


± 


0.003 


SND [15] 


-0.010 


± 


0.021+0.010 


WASA@CELSIUS [16] 


-0.026 


± 


0.010 + 0.010 


WASA@COSY [17] 


-0.027 


± 


0.008 ± 0.005 



Table 1: Theoretical predictions and experimental hndings on the slope parameter a. 

= jV n |l + az + pz z l 2 sin(30) + jz 2 + . . . } . (2.8) 

The relations between the expansion parameters up to first order in isospin breaking are 
found to be 

N c = N c x N c , N c = 1 + a{R n - R c ) , 

5 + 2Qgg - R c )b - p2 b + 3(^-^)7 J op2 rf + g(^n-^ c ) 
a = -i? c ^ , 6-i? c ^- , d = 3R c ^ , 

/= -^f' ^ = -T' «=2^' ^ = 1^' 7 = ^7- (2.9) 

The expansion in powers of R n — R c ~ 3.35 x 10~ 3 GeV 2 hinges on the fact that we have 
considered isospin breaking corrections in the definition of y (in the isospin limit, R n = R c , 
we reproduce the results derived in Ref. [6]). The relations to the Dalitz plot parameters 
of the squared value of the respective amplitudes Eq. (2.5) are then easily shown to be 

a = 2Re(a) , b = \a\ 2 + 2 Re(6) , d = 2Re(J), / = 2 Re(a6* + f) , 

g = 2Re(ad* + g) , a = Re(a) , p = Re(/9) , 7 = Re(x> . (2.10) 

The A/ = 1 rule Eq. (2.7) gives rise to relations between Dalitz plot parameters and 
normalizations of the neutral and the charged decay amplitude, namely 

N n = -3A7; , a = i(6 + 3d) . (2.11) 
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3. The modified non-relativistic effective field theory framework 

In this work we will use the modified non-relativistic effective field theory (NREFT) frame- 
work to analyze the final-state interactions in r] — > 3ir. This framework provides a useful 
tool to investigate low-energy scattering and decay processes: it has found applications in 
detailed studies of cusp effects in K — > 3ir [29-31] and r/ — ?• 37r [32] as well as r/ — > rjirir [33] 
decays, and has recently been extended to describe near-threshold pion photo- and elec- 
troproduction on the nucleon [34,35] (for an overview on cusp effects in meson decays, see 
Ref. [36]). 

An analysis of rj — > 3ir within the non-relativistic framework is useful for the following 
reasons. While the non-relativistic amplitude is perturbative, just as the chiral amplitude, 
it allows for a more accurate implementation of tttt interactions due to the inclusion of phe- 
nomenological threshold parameters as determined from Roy equations. Non-perturbative 
treatments, as for example dispersive analyses, are expected to yield yet more precise 
results. Compared to such numerically very involved studies, however, the NREFT cal- 
culation leads to a very transparent analytic representation. Moreover, it allows for the 
direct implementation of isospin breaking in particular in all kinematic effects, which is 
much more involved in ChPT and unexplored in dispersive analyses. 

In that context it is useful to narrow down the precise definition of the term "non- 
relativistic" as it is used in our work. Our representation of the decay amplitude is only non- 
relativistic in the sense that inelastic thresholds outside the physical region are subsumed 
into point-like effective coupling constants. Inside the physical region, however, we arrive 
at a fully covariant expression with the correct non-analytic low-energy behavior. The 
number of low-energy Dalitz plot couplings to be included in the Lagrangian at tree-level 
is modeled after the traditional (experimental) Dalitz plot expansion, which seems to yield 
a rather good description of the experimental data in the center of the Dalitz plot. We 
note again, see Sect. 2, that the full Dalitz plot is not accurately described by a polynomial 
expansion, since such a representation neglects non-analytic effects, such as cusps at the 
opening of the charged pion threshold (see also Refs. [8,30,32]). 

In fact, the non-relativistic approach to ij — > 3ir is not new. In Ref. [32] the authors 
performed a fit to experimental data in an attempt to investigate the cusp effect in t] — > 3ir° 
generated at the opening of the charged pion threshold. The scope of our work is entirely 
different. We focus specifically on an analysis of the Dalitz plot parameters based on 
numerical input parameters derived from ChPT. For that endeavor the amplitudes are 
calculated to yet-higher accuracy in order to ensure the incorporation of the most prominent 
effects generated by the final-state interactions. In the following section we give a brief 
introduction to the modified non-relativistic framework. 

3.1 Power counting (1): basics and tree amplitudes 

A Lagrangian treatment of r\ — > 3ir in the non-relativistic framework is provided in Ref. [30] 
and will not be repeated here. Instead, we will briefly comment on the power counting and 
outline the basic features of the amplitudes of rj — > 3tt and the mr final-state interactions. 
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A consistent power counting scheme for the modified non-relativistic effective field 
theory is constructed by introducing the formal non-relativistic parameter e and count 

• pion 3- momenta (in the rj rest frame) as 0(e), 

• kinetic energies Tj = p° — Mj as 0(e 2 ), 

• masses of the particles involved as 0(1), but A T = M 2 — M 2 as 0(e 2 ), 

• and the excess energy Q n / C = Y2i Ti as 0(e 2 ). 

Loop corrections in the perturbative series involve irir rescattering at not-too-high energies, 
which can be related to the effective range expansion of the tttt amplitude. Since these ef- 
fective range parameters are phenomenologically small, we use them as an additional power 
counting parameter, referred to generically as a wn . We thus have a correlated expansion 
in a n7r and e and can uniquely assign powers to our loop expansion (for a more detailed 
introduction to the modified non-relativistic effective field theory we refer to Refs. [29-31]). 

Following the previous counting scheme a Lagrangian framework can be constructed. 
From that the rj — > 3ir amplitude at tree level can be derived as 

■KT(si,S2, s 3 ) = K + K 1 [(p° - M^f + (p° - M n0 ) 2 + (pg - M n0 ) 2 ] + 0(e 6 ) , 
MT c (si,s 2 , s 3 ) = L + Li(p§ - + L 2 (pl - M n o) 2 + L 3 (p? - p° 2 ) 2 + 0(e 6 ) , (3.1) 

where the low-energy couplings Ki, L; L are of 0(1) and are related to the traditional Dalitz 
plot, see Sect. 3.3. The isospin relation Eq. (2.7) translates into 

K = -(310 + LiQn- L 3 Q 2 n ) , K x = -(L 2 + 3L 3 ) . (3.2) 

The number of constants included here corresponds to expanding the Dalitz plot up to 
quadratic order; we briefly comment on the possible inclusion of cubic terms at tree level 
in Sect. 4.2. We remark that the number of four independent tree-level couplings (in the 
isospin limit) chosen here equals the number of subtraction constants in several of the 
dispersive analyses [21,22] (compare Ref. [23], though). Analogously, the tttt scattering 
amplitude can be determined. We consider the following final-state processes (i) (iT a "K b — > 
ir c ir d ): (00) (00; 00), (x) (+-;00), (+0) (+0;+0), and (+-) (+-;+-). Up to 0(a^e 2 ) 
the threshold expansion of the amplitudes in the respective channels are given as 

ReT°° H = 20oo + 2A)o(s - O + ^oo(s - s$f) 2 + ^ 2 x J + .{s) + . . . , 

ReT% R = 2C X + 2D x (s - sf") + 2F x (s - sf*) 2 + ... , 

ReT+° = 20+o + 2D +0 (s - s%) + 2F +0 (s - s%) 2 - E +0 (t -«) + .. . , 

Re T+ R = 20+_ + 2D+_ (s - s^L) + 2F+_ (s - s^L) 2 -£+_(<- u) + .. . , (3.3) 

where s* hr denotes the threshold of the pertinent channel, s 1 ^ = 4M 2 , s x hr = 4M 2 , 
•s+o = (^7r° + M n ) 2 , s+I = 4M 2 . The one-loop function of the non-relativistic theory, 



i 



4M 2 



-8- 



is responsible for a cusp structure in the (OO)-channel (see Refs. [29,37] for further details). 
The low-energy couplings are matched to the effective range expansion in the following 
section. 

3.2 Matching (1): 7T7t scattering 

We want to make more sense of the low-energy couplings introduced in the previous section. 
To determine the matching relations for the low-energy constants of tttt scattering, we resort 
to the effective range expansion of the tttt scattering amplitude, which is conventionally 
decomposed into partial waves according to 

T T (s,t) = 32^(2/ + l)tf( a )Pi(z) , (3.5) 
l 

where ti(s) is the partial wave amplitude of angular momentum I and isospin /, Pi{z) 
are the Legendre polynomials, and z = cos 9 is the cosine of the scattering angle in the 
center-of-mass system. Close to threshold the partial wave amplitude can be expanded in 
terms of the center-of-mass momentum q 2 = q 2 (s) = (s — 4M 2 .)/4, leading to 

Ret? (s) = q 2l {a{ + b{q 2 + c{q A + 0(q 6 )} , (3.6) 

where a\ is the scattering length, b\ is the effective range, and cf is the (leading) shape 
parameter. In the following we use the simplified notation ai, bi, ci, as only S- and P-waves 
will be considered. In the language of NREFT power counting the previous equation is 
an expansion in orders of e, since q 2 oc e 2 . The effective range expansion is thus naturally 
related to the non-relativistic tttt scattering amplitude in Eq. (3.3), and we can read off the 
matching relations for the low-energy couplings, shown here for simplicity in the isospin 
limit: 



Coo 


16-7T 

= — (ao + 2a 2 ) , 


D 00 


4-7T 

= —(b + 2b 2 ) , 


Foo — 


|(co + 2c 2 ) , 


(3.7) 


C x 


167T 

= — (-ao + a2) , 


D x 


= y(-6o + 6 2 ) , 


F x = 


3 (-co +C2) , 




C+o 


= 8irci2 , 


D+o 


= 2vr6 2 , 


F+o = 


7T 

2 C2 ' 


E +0 = 12vrai , 


C+- 


= y(2ao + a 2 ) , 


D+- 


2-7T 

= Y (2b + b 2 ), 


F + - = 


^(2c + c 2 ) , 


E + - = 12vrai . 



Isospin-breaking corrections to these matching relations are discussed in Appendix A.l. 
Note that Eq. (3.7) is only valid up to 0(a 2 n ), i.e. tttt scattering to one loop, or r] — > 37r 
to two loops. At higher loop orders, the low-energy couplings Di and Fi are renormalized, 
which we will briefly discuss in the context of higher-loop resummation at the end of 
Sect. 4.1. 

We will use two sets of phenomenological values for the tttt effective range parameters, 
the combined Roy equation plus ChPT analysis of Refs. [38, 39] (henceforth denoted by 
ACGL) and a combination of forward dispersion relations and Roy equations [40] (KPY). 
The central or "best" values for S- and P-wave scattering lengths and effective ranges as 
obtained in those two analyses are quoted in Table 2. The determination of the shape 
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parameters is a little more delicate. We use the respective parameterizations of the phase 
shifts given in Refs. [38, 40] and calculate the scattering amplitude according to 

' 0, 2 . (3.8) 



Re^(g 2 ) = (l + -f 



1/2 tan£ 7 
1 + tan 2 Si 

Since the shape parameters are numerically very small in comparison to effective ranges and 
scattering lengths, they are rather sensitive to the method by which they are determined. 
For example, one receives rather different re- 
sults when extracting the shape parameter from 
a strict threshold expansion of the amplitude, or 
from a fit over a certain low-energy range, mini- 
mizing the x 2 -function 

X 2 ( CI ) = (W^ 2 )-g/-M 2 -c /( / 4 ) 2 , (3.9) 

in a range from the threshold 4M 2 up to the 
expansion point s n . Furthermore, the inclusion 
of an additional term diq 6 causes significant de- 
viations in the 1 = channel, since this term 
and the leading shape parameter are of compa- 
rable size. We decide to use the central values 
obtained from the minimization of Eq. (3.9) as 
the most reasonable approximation to the true 
partial wave. The numerical results for co,2 thus 
obtained are also given in Table 2. In the follow- 
ing, we use the variation between the central values of the two parameterizations [38, 40] 
as a means to estimate the uncertainty due to tttt rescattering. 

3.3 Matching (2): rj ->• 3tt 

We compare Eqs. (2.8) and (3.1) to derive the matching relation between the low-energy 
couplings of the r) — > 3ir tree amplitude and the traditional Dalitz plot parameterization, 
namely 





ACGL 


KPY 


a 


0.220 


0.223 


a 2 


-0.0444 


-0.0444 


b xM% 


0.276 


0.290 


b 2 x Ml 


-0.0803 


-0.081 


c x 10 2 M 4 


-0.19 


0.04 


c 2 x 10 2 M 4 


1.33 


0.68 


ai x 10M 2 


0.379 


0.381 


6i x 10 2 M 4 


0.567 


0.512 



Table 2: Input values for the scatter- 
ing lengths aj, effective ranges bj, and 
shape parameters cj as determined from 
the two parameterizations ACGL [38, 39] 
and KPY [40] (see text for discussion). 



L = AA c trcc (l + ~a tTCC R n + b tTCC R 2 n ) , 
L 2 = 4AA c trce M 2 6 tree , 



K x = 4A^ ree M 2 a tree , 

L\ = -2Nl rcc M v {a trce + 2b trcc R n ) , 

L 3 = 4AA c trcc M 2 ci tree , (3.10) 



where the superscript "tree" denotes tree-level input parameters. Note that Eq. (3.10) 
fulfills the isospin relation Eq. (3.2) as long as Eq. (2.11) is satisfied. To extract the Dalitz 
plot parameters in the non-relativistic framework, we have to fix the numerical input for 
the tree-level low-energy couplings for the r] — > 3tt amplitude. We determine the low-energy 
couplings of the Dalitz plot in Eq. (3.10) by matching the non-relativistic framework to 
the one-loop ChPT amplitude [5] at the center of the Dalitz plot. Following Ref. [8], we 
evaluate the chiral rj — > 3tt amplitude using neutral masses everywhere. 
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We remark that the upcoming [22] (and previous [21]) dispersive analyses use the Adler 
zero of the rj — > 7r + 7r _ 7r° amplitude as the matching point, compare Eq. (2.6). It is protected 
by SU(2) symmetry and therefore not prone to large strange-quark-mass corrections. The 
chiral series is thus expected to converge rather quickly, which makes the Adler zero a 
natural choice. The fact that it lies outside the physical region at roughly sa ~ ^M 2 , 
however, renders matching the non-relativistic framework to the chiral amplitude at this 
point ill-fated: the expansion in terms of e does not necessarily converge there, and we 
therefore have to resort to matching inside the Dalitz plot. 

For the matching procedure we tune the rescattering parameters in the non-relativistic 
amplitude in such a way as to mimic the chiral amplitude. In essence this means that 
the scattering lengths and effective ranges are fixed at their current algebra values (this 
corresponds to the insertion of 0(p 2 ) vertices in the chiral expansion). Explicitly, we have 

7 'Ml n A Ml 



CA _ it CA _ 7v CA 



32vrF2 ' 2 16vrF2 ' 1 24vrF2 ' 

ACA _ 1 h CA l _ fo nl 

We proceed analogously with the 77 — > 37r couplings that enter the non-relativistic amplitude 
at one-loop level and derive from Eq. (2.6) 

(M 2 -M 2 )(M 2 + 3M 2 ) 3 
c " 16QV3X 2 M2 ' ~M 2 -M 2 ' [6AZ) 

For our numerical analysis we will use the value for Q dictated by Dashen's theorem, 
Qd = 24.2. Note that the specific choice does not hold any ramifications for our main 
statements, since it merely enters in the normalization, which drops out in the Dalitz plot 
parameters. 

The above matching procedure is consistent as it ensures that the imaginary parts 
are exclusively generated by tttt final-state interactions. Residual effects from the chiral 
pion loops are purely real and absorbed in the low-energy couplings. We use matching to 
C(p 4 ) and not to 0(p 6 ) for practical reasons: the above matching procedure is simpler and 
our results can be used to compare with and interpret the dispersive analyses directly. A 
high-precision determination of the Dalitz plot parameters would likely require matching 
to 0(p e ), but for that purpose the low-energy constants showing up at 0(p & ) may not 
be known with sufficient accuracy. Numerically we obtain from matching to the ChPT 
amplitude at 0(p 4 ) (using the chiral SU(3) low-energy constant 2 L\ = -3.5 x 1CT 3 [41]) 

A? c trec = -0.158 , a trcc = 13.428 GeV" 2 , 

b trcc = -7.291 GeV" 4 , d trcc = 5.189 GeV" 4 . (3.13) 

The particle masses used throughout this work are given by the current particle data 
group values [42], i.e. M n = 139.57 MeV, M w o = 134.98 MeV, and M v = 547.86 MeV. 



2 The effects of varying L\ within its error were checked to be tiny compared to other uncertainties. We 
therefore only use the central value. 
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3.4 Power counting (2): loops and rj — > 3ir 

The power counting scheme discussed in Sect. 3.1 gives rise to a natural decomposition 
of the NREFT amplitude. This can be seen as follows. The modified non-relativistic 
propagator counts as 0{e~ 2 ) (see e.g. Ref. [29]), the loop integration measure (with one 
energy and three momentum integration variables) as e 5 , therefore any loop integral with 
two-body rescattering contributes at 0(e). Moreover, such a loop always involves a tttt 
rescattering vertex and is thus of 0(a n7r ). The decomposition of the full r\ — > 3-ir amplitude 
according to its loop-structure, 

M n/c ( Sl , 8 2 , S 3 ) = MllJ C c ( Sll 8 2 , S 3 ) + M 1 J°° P (S 1 ,82, S 3 ) + M^ P '(«!, 8 2 , 8 3 ) + . . . , (3.14) 

is thus an expansion in powers of a n7T e. There is an interesting simplification of Eq. (3.14) 
close to the center of the Dalitz plot (s\ « S2, S3 ~ s n ) above all two-pion thresholds. The 
contribution of the one- loop function is purely imaginary as can be seen from Eq. (3.4). At 
the same time the two-loop bubble diagram, which is the product of two one-loop functions, 
is purely real and it can be shown that the imaginary part of the non-trivial two-loop 
function does not contribute at this order (see for example Refs. [29,30]). Symbolically we 
can write both amplitudes in terms of the power counting parameter a nw , 

M = TWtrcc + i-Mi-ioopa™ + 7W 2 -ioo P oL + 0(zo^e 3 , ia^e 4 ) , (3.15) 

where the ©(ia^e 4 ) term stems from the three-particle cut at two-loop order, which is 
numerically small as discussed in Appendix C and therefore neglected. By taking the 
absolute value squared we obtain 

\M\ 2 = M 2 tICC + (M 2 _ loop + M tICC x M 2 . loop )a 2 n7T + 0(a^e 4 , a 3 n7T e 5 ) . (3.16) 

We therefore expect one- and two-loop effects to be of the same size at the center of the 
Dalitz plot, as only the two- loop contributions can interfere with the dominant tree terms 
there, and thus to impact the Dalitz plot parameters about equally. 

The heightened importance of rescattering effects in Dalitz plot parameters is further 
substantiated by another observation. Consider the generic one-loop function of tttt rescat- 
tering in the non-relativistic theory expanded about the center of the Dalitz plot (s = s n , 
we neglect isospin-breaking effects in the following discussion, so that M n o = M n ): 



T , , _ A 6Mg s-s n 18M 2 (s n -3M 2 ) / s - s n \2 

(S) 16vr V + 8 „ M2-9M2 si \Ml-^Ml) 1 ' 



« 77 u±y± ix °n ±y± r) 

= 0(e) , (3.17) 

since s — s n = 0(e 2 ) and M r] — 3M W = 0(e 2 ). The same holds true for the two- loop 
functions. This implies that contributions to higher-order Dalitz plot parameters from the 
loop functions are enhanced non-analytically in M v — 3M n . We conclude from Eq. (3.16) 

•^-loop + -Mtrcc x M 2Aoop = 0(e 2 ) , (3.18) 
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Figure 2: The Feynman graph topologies at one and two loops contributing to the decay rj — > 'Sir 
in NREFT. The double line denotes the r) particle, the single lines stand for pions (of arbitrary 
charges) . 

which has substantial consequences for the slope parameter a of the neutral decay channel. 
We can parameterize the slope parameter according to 



From relations Eq. (3.10) we find d tree = 0(1) and consequently the slope parameter 
at tree-level is of order ao oc Q 2 d trcc = 0(e 4 ), whereas rescattering effects enter the slope 
parameter at 0(a 2 r7r e 2 ). This obviously implies that rescattering effects become increasingly 
more important for higher-order Dalitz plot parameters. On the other hand, they are 
far less significant (as we will confirm numerically below) for the normalization of the 
amplitude, for which we expect higher-order quark-mass renormalization effects to be more 
important. 

The full NREFT representation beyond tree level with isospin breaking included is 
given in Appendix B. It comprises the loop graph topologies displayed in Fig. 2, and is 
fully consistent in terms of non-relativistic power counting up-to-and-including 0(a 2 n e 4: ), 
i.e. the vertices of the two- loop graphs are included at 0(e 2 ). Phenomenologically, one finds 
that the expansions of the r\ — > 3tt and inr — > irir polynomials in powers of e 2 only converge 
well starting from next-to-next-to- leading order, i.e. the 0(e 2 ) terms (the linear slope in 
rj — > 7T + 7r~7r and the -kit effective ranges) are not really suppressed compared to the leading 
(constant) terms. This observation is readily understood resorting to chiral perturbation 
theory: due to the Goldstone nature of the pions, the constant terms are chirally suppressed 
by powers of M 2 and the leading 0(p 2 ) amplitudes are linear in energy s. In other words, 
the 0(e 2 ) contributions are "suppressed" versus the constant ones by factors of s/M 2 , and 
only starting from 0(e 4 ), the relative suppression is s/A 2 with A x ~ 1 GeV. The chiral 
two-loop or 0(p & ) calculation [6] contains all the leading 0{p 2 ) vertices and therefore the 
linear 77 — > 37r slope as well as effective ranges for the ttit interaction (although not quite the 
phenomenologically accurate ones). In order to guarantee that our NREFT representation 
of the decay amplitude is at least as accurate as the chiral two-loop one, we include all 
combinations of linear energy dependences in the three vertices of the two-loop diagrams. 
Thus, our amplitude also contains terms that are of 0(a 2 n e ti ) and 0(a 2 n e 8 ), and due 
to the enhancement discussed above, the numerically most important ones appearing at 
those orders. The representation of the "double bubbles" (see also Fig. 2) is even strictly 
complete up to 0(a 2 7T e d ), as P-wave contributions only start at 0(a 2 7r e 8 ). Furthermore, 
we have added shape parameter terms in the "double bubbles" and in the outer vertex of 
the irreducible two-loop graph, where the addition of these terms is trivial. 



(3.19) 
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4. The isospin limit 



We first give an analytic and numerical treatment of the amplitude in the isospin limit, 
which we define as M n o = M w and using Eq. (3.2). This already includes the gross features 
of our total analysis. However, in the isospin limit, we can give relatively simple closed 
analytic expressions for all parts of the amplitudes up to two loops. 

4.1 Structure of the amplitude 

The non-relativistic decay amplitude (for the charged channel) can be split into parts 
consisting of tree and final-state contributions 

M C (S!, S 2 , S 3 ) = Mc CC (s 1 ,S 2 , S 3 ) + Mf(s!, S 2 , S 3 ) , (4.1) 

where the tree amplitude is given by 

MT( Sl , 8 2 , s 3 ) = AA c trcc {l + a trcc ( S3 - s n ) + 6 trec (s 3 - s n ) 2 + J tree ( Sl - s 2 ) 2 } , (4.2) 

and the rescattering contributions of the amplitude can be decomposed up to 0(p 8 ) ac- 
cording to the isospin structure of the final-state pions [21,43,44], 



Mf(s 1 ,s 2 , s 3 ) = M (s 3 ) + (s 3 - S!)Mi(s 2 ) + (s 3 - s 2 )M 1 (s 1 ) 

2 

+M 2 {s 1 ) + M 2 {s 2 ) - -M 2 {s 3 ) , 



(4.3) 



where the index / = 0, 1, 2 of the function Mj(si) denotes the total isospin of the respective 
kinematic channel. At 0{a^ qT e A ) (for details see Sect. 3.4) the isospin amplitudes are given 
as 



+ 



5 
3 
32tt 



M„(.s) - -\e (s)J(s)(l + 16TTa (s)J(s)) 



21 



1 q 2 Ci( s ) T , x 807rs 



+ 



(e' (s)a (s) + 2£' 2 (s)a 2 (s))F^(s) + ^(-a (5) - a 2 (s)) - i' Q {s)b 

^ s)h ) W p{1) (s) - 2Li (f 60 - b2 ) ^ F(2) {s) \ } i67ra ° (s) ' 

(e' (s)a (s) - £' 2 (s)a 2 (s)) (F^(s) - 2F^(s)) 

^(|°o(S) + - £' (s)b + t' 2 {s)b 2 ) ^V (1) 00 - 2F^(s)) 

M„Q 4 



2Q C 
4irai(s) 



M 2 (s) = U 2 (s)J{s)(l + 16ira 2 (s)J(s)) 
16vr 



%(s)a (s) +e' 2 (s)a 2 (~s))F(°\s) + [^(a (~s) - ^) - 5t' (s)b 
-^(^ 2 )^F« (s) -4L 1 (6 -^)^F( 2 )( S )]}l6 7 ra 2 ( S ) , (4.4) 
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where the various polynomials are given by 



4>(s) 


= -J(s) + -Ms) 




£ 2 {s) 




h(s) 


= L 1 + 2L 2 (^-- 


£'{s) 


- Lo + Ll ( 2 QO- 


ai(s) 


= a/ + + ciq 



(s) = L + Lx (p° - M,,.) + L 2 (p° - M,) 2 + L : 



Q 2 



U + 2L 2 (^ - M n ) + 2L 3 (p° - ^) , e' (s) = + , 



+ ^3 



4Q 2 . 

2 A2 



, Q 2 



'M 

I + (-^ M„ 



7 4 , 

and we use the kinematic variables 

Q _M* + MZ-s M 2 -M 2 + s 



(4.5) 



2M„ 



2M„ 



X(M 2 ,M 2 ,s) 



s ~ = 2M 2 - S + -^(. + 2Q 2 ) , 



(4.6) 



with the Kallen function X(x,y,z) = x 2 + y 2 + z 2 — 2(xy + yz + zx). Note that the 
shape parameter terms oc cj are to be omitted in aj(s); we also neglect them in the 
1=1 partial wave. We use the shorthand expressions J(s) = J + -(s) and F( n \s) = 
F^ 1 \m w , M n , M n , M n , s) (in the isospin limit), where for the exact form of the two-loop 
functions we refer to Appendix B. We can now write the Dalitz plot parameters in terms 
of the isospin amplitudes, namely for the charged channel 



M c = Mr c + Mo(s n ) + -M 2 (s, 



a = - j| (A? c tree « tree + M$\s n ) + 3Mi(a n ) - j^W) , 



J\ c \ 



i ± yo nJ -r ^/v l2 



8AA C 

and for the neutral channel 



« = -^{Mf\s n )-Mf (s n )), 



a 



4A/" C 



(2A? c tree (fe tree + 3d trcc ) + M { *\s n ) + ->4 2) (sn)) , 



(3), 



(4.7) 



(4.8) 
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Figure 3: Diagrammatic expression of the resummed amplitudes. Above: the bubble chain. Below: 
resummed external vertex of the non-trivial two-loop graph. The line style is as in Fig. 2. 

where M.\ (s n ) denotes the ra-th derivative of the function Aii(s), evaluated at the center 
of the Dalitz plot. Note that d and g do not receive contributions from the isospin 1 = 
amplitude. 

Despite working in the limit of evaluating all amplitudes for the charged pion mass, we 
employ the physical values for R c and R n in Eqs. (4.7), (4.8). These prefactors stem from 
the conversion of a, a etc. into a, a etc., see Eq. (2.9), and are just due to a normalization 
choice in the definition of the Dalitz plot variables x and y; we therefore decide to present 
our results including this "trivial" isospin-breaking effect already at this stage. Note that 
due to (M v - 3M 7r o)/(M, - 3M n ) rj 1.11 and (M v - 2M W - M 7r o)/(M ?? - 3M V ) ps 1.04, the 
effects of using these normalization factors in the isospin limit are large, most so for the 
neutral channel, where a for instance is affected by a shift of 22%. 

In our numerical analysis we will observe that among the two-loop contributions those 
of the non-trivial two- loop graphs, see Fig. 2 (right), are in general strongly suppressed. 
This can be traced back to the isospin properties of these pieces: for those Dalitz plot 
parameters to which the 1 = partial wave can contribute, it usually dominates. For 
those graphs that only describe rescattering in one channel and can be written as simple 
products of one-loop functions, see Fig. 2 (middle), the 1 = isospin amplitude receives 
contributions proportional to second powers of ao, bo, etc., whereas the "inner" vertex in 
the non-trivial two-loop contributions has parts of / = and 1 = 2 (P-waves vanish due 
to symmetry reasons in the isospin limit) that tend to partially cancel each other. In an 
attempt to estimate (partial) higher-order corrections, we therefore expect to find a good 
approximation to the full result by iterating the bubble diagrams and the exterior two- 
particle rescattering of the non-trivial two-loop function as depicted diagrammatically in 
Fig. 3. In the aforementioned representation the unitarized amplitudes are easily deter- 
mined to be 

m», ,_Mo(s)-IUs){l67ra ( S )J(s)) 2 

Jvl (S) — — 



1 — 167TOo(s) J(s) 

Mi(s) 



1 - lGir ai (s)q 2 J(s) 



1 — WTTa2{S)J{S) 

The inclusion of iterated diagrams requires modified matching relations for the effective 
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range parameters. This becomes obvious when considering the expansion of the iterated 
bubble sum of tttt scattering of isospin / = 0, 2 at the tttt threshold: 

,,3 \ , „3 _l „5 _ Q^2 A ^ /f 2, 

Re 



arts) 



1 — 16nai(s) J(s) 



= ai + ~M) q + + M4 K + } • 

(4.10) 



One immediately sees that the effective range picks up a contribution from two-loop dia- 
grams, the shape parameter from two- and four-loop diagrams. To account for this shift, 
the above expression has to be compared with the effective range expansion of the tttt 
amplitude (for 1 = 0), 

Ret I (q 2 )=a I + b I q 2 +c I q' l + O(q 6 ) , 1 = 0,2, (4.11) 

from which one reads off the following renormalization prescriptions: 

„ren _ , ren _ . , a J ren _ „ gf ~ 2a J ~ 3af6jM^ 

Oj - a/ , b 7 - 6/ + — f , C/ - c/ — 3 . (4.12) 

The af n , b r f n , c r f n are now to be inserted into the matching relations for the coupling 
constants Cj, Dj, i^. The renormalization prescriptions have pretty remarkable effects 
in the isospin 1 = channel, where the shape parameter is shifted from — 0.002M" 4 to 
-l-O.OSOM^ 4 (for the ACGL parameter set). We note that the P-wave effective range b± 
does not pick up an additional contribution due to the q 2 (s) prefactor. Corrections in the 
P-wave channel start at 0(e 8 ), that is the higher-order shape parameter d\. 

4.2 Numerical results 

We begin our numerical analysis of the various 77 — > 37r Dalitz plot parameters by investi- 
gating how the tree-level values are modified at one- and two-loop order, and finally beyond 
two loops (via the estimate through the unitarized amplitudes in Eq. (4.9)). This part of 
the analysis is based solely on the ACGL parameters for the tttt final-state interaction; the 
qualitative conclusions are identical for the KPY parameterization. We keep the rj — > 3tt 
tree level parameters fixed as obtained by matching to ChPT at 0(p 4 ) throughout, see 
Sect. 3.3. Our results are summarized in Table 3. In particular, we observe the following: 

1. Individual loop corrections to the Dalitz plot parameters are sizeable; their rela- 
tive importance grows with increasing order (in e) of the parameters concerned, as 
suggested a priori by power-counting arguments (see Sect. 3.4). 

2. One- and two-loop contributions are in general of the same size, as indeed expected, 
with a tendency to cancel to varying extent due to contributions of opposite sign. 
This again substantiates the power-counting arguments of the NREFT framework, 
which is particularly interesting in the case of a: while at one loop we see a sizeable 
positive shift added to the already positive tree-level result, the two-loop correction 
overwhelms both, leading to a negative total. We therefore find the correct sign for 
a, as opposed to the ChPT result. At two loops our result is in fairly good agreement 
with the dispersive one from Rcf. [20] . 
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full two-loop 


0.0287 
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-0.045 


unitarized 


0.0284 


-1.268 
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0.053 
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-0.042 
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tree 


0.279 


0.0107 





0.0001 






one-loop 


0.304 


0.0227 


0.0005 


0.0000 






two-loop* 


0.260 


-0.0209 


-0.0027 


0.0007 






full two-loop 


0.258 


-0.0192 


-0.0036 


0.0009 






unitarized 


0.255 


-0.0249 


-0.0043 


0.0013 







Table 3: Results for the charged and neutral Dalitz plot parameters in the isospin limit. We 
show tree, one-loop, two-loop neglecting the irreducible two- loop graphs (marked "two-loop*"), and 
full two- loop results calculated as described in the text, plus the result employing the unitarized 
amplitudes according to Eq. (4.9). 

3. There are large contributions from derivative couplings at two-loop order. This is seen 
when considering the amplitude expanded only up to 0(a W7r e 5 , a^ w e 2 ) (cf. Ref. [30]), 
at which order only constant vertices are implemented at two loops. In this approxi- 
mation, we find numerically e.g. a = +0.033. Once the effective range corrections in 
the 1 = two-loop bubble are added, a receives a shift to —0.017. This observation 
explains why the authors of Ref. [32] obtain a positive sign for a when matching to 
ChPT at tree-level: no derivative couplings at two-loop level are included in that 
work. With respect to this omission, matching to ChPT at tree-level plays a minor 
role in the deviation from our result. 

4. By comparing the two-loop contributions with and without the parts due to the 
irreducible two-loop graphs, see Fig. 2 (right), we see that at least in those parameters 
that receive contributions from the 1 = amplitude the irreducible two-loop graphs 
only give a very small contribution. As detailed before, this can be traced back to 
the isospin structure of the different amplitudes. Specifically a is a case in point: 
the simple "bubble sum" type two-loop graphs shift it by about —0.044, while the 
irreducible graphs only add +0.002. 

5. Our estimate of higher-order effects via simple two-channel unitarization shows that 
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ACGL 2-loop 


ACGL unit. 


KPY 2-loop 


KPY unit. 


average 




0.0287 


0.0284 


0.0285 


0.0282 


0.0284 ± 0.0002 


a 


-1.290 


-1.268 


-1.291 


-1.267 


-1.279 ± 0.012 


b 


0.379 


0.342 


0.382 


0.340 


0.361 ± 0.021 


d 


0.056 


0.053 


0.052 


0.050 


0.053 ± 0.003 


f 


0.071 


0.101 


0.073 


0.107 


0.089 ±0.018 


9 


-0.045 


-0.042 


-0.043 


-0.041 


-0.043 ± 0.002 


a 


-0.0192 


-0.0249 


-0.0227 


-0.0291 


-0.0242 ± 0.0049 


P 


-0.0036 


-0.0043 


-0.0043 


-0.0051 


-0.0043 ± 0.0007 


7 


0.0009 


0.0013 


0.0013 


0.0017 


0.0013 ± 0.0004 



Table 4: Results for charged and neutral Dalitz plot parameters in the isospin limit with different 
input on 7T7T scattering parameters from Refs. [38] and [40]; see Sect. 3.2. Shown are the results 
both for two loops and for the unitarized amplitudes. \Af n \ 2 = 9|7V C | 2 is not shown separately. 

those are significantly smaller than the (individual) one- and two-loop effects, al- 
though not negligible throughout. Due to the smallness of the irreducible two- loop 
graphs, we expect to catch the major part of the higher-order corrections in this way. 

In order to study the dependence of our results on the precise input for tttt scatter- 
ing, we next compare the values obtained for the various charged and neutral Dalitz plot 
parameters, at two loops and unitarized, for the ACGL and the KPY parameter sets in 
Table 4. In most cases, the variation with different irir input is a bit smaller than the 
difference due to the higher-order estimates, although not by much. As our final result in 
the last column of Table 4, we determine central values and (symmetric) errors in such a 
way as to cover all four values for each parameter. 

Our finding for the r\ — > 37r° slope parameter, a = — 0.024±0.005, is considerably closer 
to the current experimental average a = —0.0317 ± 0.0016 [42] than previous theoretical 
approaches. Note again that the theoretical prediction for a is lowered (in absolute value) 
by about 22% if the charged pion mass is used in the definition of z. We predict the 
(yet unmeasured) higher-order Dalitz plot parameters /3 and 7 in the neutral channel 
to be different from zero, but very small. In particular, neglecting a term oc 7Z 2 in an 
experimental extraction of a based on the radial distribution dT/dz alone (in which a term 
oc f3 cancels for z < 0.756, compare Fig. 1) should affect a by less than the value of 7, 
hence still below the current uncertainty, although not by much given the precision of the 
most recent experimental determinations. 

As we will see below, there are sizeable isospin-breaking shifts in the charged Dalitz plot 
parameters. We therefore defer a detailed comparison to experimental values to Sect. 5. 
We only wish to make a remark on the cubic parameters / and g here. Apart from the 
fact that a large contribution to these is given by 2Re(a6*) and 2Re(ad*), respectively, the 
remainders (or /, g) are given entirely in terms of loop contributions. If we, in addition, 



- 19 - 



allow for cubic tree level terms f tiee and g tree and match the latter to ChPT at 0(p 4 ), the 
total results receive shifts of —0.002 and —0.011 hence very and relatively small effects, 
respectively. Although chiral 0(p 6 ) corrections might modify these numbers significantly, 
we still regard them as indications that the dominance of loop contributions (as suggested 
by e power counting) holds here. 

4.3 Comparison to a in ChPT at two loops 

While dispersive analyses find values for a similar to ours [20], a serious puzzle is the 
question why the calculation of this quantity in ChPT to two loops [6] does not arrive 
at least at a negative value for a "naturally", i.e. as the central value (disregarding the 
large error bar due to the estimated fit uncertainty). After all, in addition to potentially 
significant chiral SU(3) renormalization effects of what would be subsumed in the tree-level 
couplings of the NREFT representation, ChPT at 0(p 6 ) also includes all the pion two-loop 
graphs shown to be important here. 

It turns out that this failure of the chiral two-loop calculation can partly be understood 
within our framework, investigating rescattering effects only, but of course neglecting the 
0(p 6 ) modified tree-level couplings. In order to mimic the chiral expansion, we note that 
in an 0(p 6 ) calculation, the tttt vertices inside two-loop graphs are only included to their 
current-algebra (or 0(p 2 )) accuracy, see Eq. (3.11), while inside the one-loop diagrams, tttt 
rescattering is taken care of up to 0(p 4 ). By inserting the respective values for the tttt 
threshold parameters in our amplitude, we find 

a ChPT = -0.0011 , (4.13) 

hence a value close to zero. We attribute the remaining difference to the central result for 
a in Ref. [6] to different tree-level couplings as determined in that paper. As we found 
that precisely the two-loop effects turn a negative, dominated by the 1 = amplitude, we 
conclude that a large part of the discrepancy between ChPT at 0(p 6 ) and our result (or the 
one from dispersion relations) is due to the significantly weaker tttt rescattering (compare 
e.g. clq A ~ 0.16 vs. ao = 0.220 from Ref. [38], which enters the two-loop effects squared). 
The precise choice of the set of rescattering parameters therefore has a large effect on the 
result for a (and, slightly less dramatically so, on other Dalitz plot parameters). The 
inclusion of improved values for the effective ranges and shape parameters produces a large 
shift of the chiral result towards the experimental value. 

In a very condensed manner, we can therefore point to one specific diagram, Fig. 2 
(middle), which accounts for roughly half of the discrepancy between the central value of 
the chiral prediction at C(p 6 ) and the experimental value for a. More specifically, the 
discrepancy is caused by contributions of the diagrammatic topology of this kind. Since 
(at least) next-to-leading order contributions to the tttt vertices are required, one needs to 
include these diagrams up to 0(p 8 ) and higher in strict chiral power counting. To substan- 
tiate this claim and ensure that it is not an artifact of the non-relativistic framework, we 
replace the non-relativistic two-point function J(s), Eq. (3.4), by its relativistic counter- 
part J- K - K {s), which differs from the former by its real part (given explicitly in Eq. (A. 11)). 
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Doing so requires a different matching procedure to account for the (otherwise absent) mass 
renormalization effects on the various coupling constants thus induced; we will not spell out 
this exercise in detail. The main conclusion however is fully consistent with our findings 
above: the "double bubble" graphs alone shift a by —0.042 (to be compared with —0.044, 
see Table 3); calculating them with current algebra values for the tttt threshold parameters 
reduces this effect by nearly a factor of two, which corresponds to the discrepancy between 
Eq. (4.13) and the value obtained in NREFT. 

One might argue that a parameter as subtle as a could also be subject to other very 
sizeable 0(p e ) corrections; in particular, contributions from chiral low-energy constants 
appear for the first time at that order. For a superficial impression of these effects, we 
investigate precisely the 0(p®) polynomial in the amplitude calculated in Ref. [6]. One 
easily finds the following combination of low-energy constants contributing to a: 



a 



(6) 
LEC 



2 



F 4 

1 TT 



± {CI + CI + 3C r 9 + CT - 2C[ 2 + 2C r 22 + 3C r 24 + C r 25 ) . (4.14) 



The couplings Cf are estimated in Ref. [6] using resonance saturation. Vector contributions 
cancel in Eq. (4.14), as they must, with no P-waves appearing in the neutral decay channel. 
Using the scalar resonance estimates given in Ref. [6], we arrive at the very simple and 
compact expression 

4 6 ec = « 0.005 , (4.15) 

where c m = 0.042 GeV, c d = 0.032 GeV, and M s = 0.98 GeV. There are serious doubts 
about the reliability of the resonance saturation hypothesis in the scalar sector [45]; indeed 
on might argue that the masses of even heavier scalar states ought to be used in Eq. (4.15), 
further suppressing their contribution to a. We nevertheless confirm that contributions 
from the ChPT low-energy polynomial at 0(p 6 ) are rather small; in particular they have 
a positive sign, so they cannot serve as an alternative explanation to arrive at a negative 
a. We also emphasize that the above is only a very rough estimate of the expected size of 
the effects and does not by any means replace a consistent matching procedure. 



5. Isospin breaking in r] — >• 3ir 

In this section we discuss higher-order isospin-breaking contributions to the decay rj — ¥ 3ir. 
We concentrate on the following four contributions: 

1. Isospin breaking in i] — > 7r + 7r~7r° due to Q n / Q c . There are significant corrections to 
the charged Dalitz plot parameters due to the terms oc (R n — R c ) in Eq. (2.9), which 
stem from the subtleties in the definition of the center of the Dalitz plot discussed in 
Sect. 2. 

2. Other isospin corrections due to the difference between the charged and the neutral 
pion mass. These in particular concern the incorporation of the correct thresholds 
inside the loop contributions, which is necessary for a description of the boundary 
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regions of the Dalitz plot, among them the cusp effect in r/ — > 3ir . The represen- 
tation of the amplitude in the non-relativistic framework allows us to work in the 
particle (and not in the isospin) basis, and thus we can incorporate mass effects in a 
straightforward fashion. 

3. Isospin-breaking corrections to the tttt rescattering parameters. We use the phe- 
nomenological values for the scattering lengths and effective ranges, which have been 
determined in the isospin limit [38,40], and calculate corrections to each channel 
from the one-loop SU(2) 7T7r-scattering amplitudes with electromagnetic corrections 
included. 

4. Next-to-leading-order isospin-breaking effects in the 77 — > 37r tree level couplings, 
calculated in one-loop ChPT, which modify Eqs. (2.11). 

The representation of the NREFT amplitude to two loops with fully general masses and cou- 
pling constants, allowing for all of these isospin-breaking effects, is given in Appendix B.l. 
Furthermore Appendix B.2 shows the generalization of the unitarization prescription given 
in Eq. (4.9) for the case of isospin violation. In our numerical evaluation we will add these 
contributions cumulatively to the results of Sect. 4. 

In this context we should comment on radiative (real- and virtual-photon) correc- 
tions to these decays. In order to be able to sensibly discuss a Dalitz plot expansion 
of the squared amplitudes in question, we assume that the universal radiative corrections 
(Gamow-Sommerfeld factor, bremsstrahlung contributions etc.), as discussed in the frame- 
work of NREFT in Ref. [31], have already been subtracted from the experimental data when 
determining Dalitz plot parameters. In order to extract the corrections of point 4 above 
from the calculation in Ref. [8], these subtracted contributions have to be matched cor- 
rectly, as detailed in Appendix A. 2. The non-universal or "internal" radiative corrections 
that play an important role in the analysis of the cusp effect in K — > 3ir [31,46] do not 
have a similarly enhanced effect in the center of the Dalitz plot. From the point of view of 
chiral power counting of isospin-breaking corrections, these constitute higher-order effects 
than those considered consistently in Ref. [8] (as they only appear at two loops); further- 
more, in 7] — > 7r + 7T~7r°, even diagrams beyond those calculated in Ref. [31] would have 
to be included. We have checked, though, that the effect of photon exchange inside the 
charged-pion loops on the rj — > 3ir° Dalitz plot expansion is small, even on the scale of the 
other small isospin-breaking effects discussed below. 

The by far largest isospin-breaking effects on the Dalitz plot parameters, beyond the use 
of the correct overall normalization factors of Q n and Q c in the definitions of the kinematic 
variables that was already incorporated in the previous sections, are the kinematic effects 
due to the fact that for the decay r] — > it + it~it the position defined by x = y = does 
not coincide with s\ = S2 = S3 when M n / M n o . Using the correct prescriptions given in 
Eq. (2.9), we find the results listed in the left column of Table 5. The corrections are very 
sizeable: our analysis shows that a is reduced (in magnitude) by 5%, b by even 14%. These 
kinematic effects constitute the bulk of the isospin breaking corrections to the charged 
parameters. 
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0.0310 ±0.0003 


0.0309 


± 0.0003 


0.0310 ± 0.0003 


0.0310 ± 0.0003 


a 


-1.218 ±0.013 


-1.214 


± 0.013 


-1.214 ± 0.014 


-1.213 ± 0.014 


b 


0.314 ±0.023 


0.310 


± 0.023 


0.308 ± 0.023 


0.308 ± 0.023 


d 


0.051 ±0.003 


0.051 


± 0.003 


0.050 ± 0.003 


0.050 ± 0.003 


f 


0.084 ±0.019 


0.082 


± 0.018 


0.083 ± 0.019 


0.083 ± 0.019 


9 


— U.Uoy ± U.UUz 


—U.Uoy 


± U.UUz 


n nQn _i_ n nno 
— U.Uoy ± U.UUz 


n nQn _i_ n nno 
— U.Uoy ± U.UUz 


\M n \ 2 


0.256 ±0.002 


0.256 


± 0.002 


0.255 ± 0.002 


0.252 ± 0.008 


a 


-0.0242 ± 0.0049 


-0.0241 


± 0.0049 


-0.0247 ± 0.0048 


-0.0246 ± 0.0049 


P 


-0.0043 ±0.0007 


-0.0043 


± 0.0008 


-0.0042 ± 0.0007 


-0.0042 ± 0.0007 


7 


0.0013 ±0.0004 


0.0012 


± 0.0004 


0.0013 ± 0.0004 


0.0013 ± 0.0004 



Table 5: Central results for the charged and neutral Dalitz plot parameters with isospin breaking 
in kinematic relations, masses in loop functions, irir threshold parameters, and 77 — > 3ir° tree level 
couplings (see text for details). 



The modifications that arise from using physical pion masses in the loop functions and 
derivative couplings are very small in the expansion around the center of the Dalitz plot, 
see Table 5 (second column). The charged parameters are typically reduced in magnitude 
on the level of about 1%; a is shifted by ±0.0001 only, an order of magnitude below the 
uncertainty due to different irir parametrizations. The importance of pion-mass effects in 
loops only becomes really visible when studying the full Dalitz plot distribution also at its 
boundaries (see Refs. [8,32]). 

The next column in Table 5 shows the effect of isospin-breaking corrections in the 
7T7T threshold parameters. For this purpose, we have calculated the electromagnetic con- 
tributions to the matching relations up-to-and-including 0(e 2 p 2 ) in the chiral expansion 
for S- and P-wave scattering lengths and S-wave effective ranges, using the results for the 
one-loop 7T7T scattering amplitudes in the presence of virtual photons of Refs. [47,48]. The 
necessary matching procedure is described in detail in Appendix A.l. The modifications 
in the Dalitz plot parameters are largest for a, where a 3% effect is observed. The contri- 
butions to the remaining parameters stay well below or around 1%. In all cases the shifts 
are dominated by isospin-breaking corrections in the S-wave scattering lengths and thereby 
the 0{e 2 ) chiral corrections, as expected by power counting. 

Finally, we want to investigate the effects of isospin breaking on the relations in 
Eq. (2.11), i.e. next-to-leading order isospin breaking in the 77 — > 3ir tree level couplings. 
These can be extracted from the chiral one-loop calculation of the 77 — > 3tt decay amplitudes 
to 0(e 2 (md — m u )) in Ref. [8]. We write the corrections in the form 

Af n = -3tf c + A#, a=i(6 + 3d) + A fi , (5.1) 

where = 0(e 2 (m ( i — m u )) and A5 = 0{e 2 ). Note that no corrections of 0((m<2 — m u ) 2 ) 
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Crystal Barrel [49] 
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± 


0.07 


0.22 


± 


0.11 


0.06 


db 


0.04 (input) 


Layter et al. [50] 


-1.08 


± 


0.014 


0.034 


± 


0.027 


0.046 


± 


0.031 


Gormley et al. [50] 


-1.17 


± 


0.02 


0.21 


± 


0.03 


0.06 


± 


0.04 



Table 6: Results for the charged Dalitz plot parameters in comparison with various theoretical 
and experimental determinations (the next- to- leading-order errors are only due to L\). 



(in A_^) and 0(md — m u ) (in A„) occur, respectively. The analytic results of the expansion 
and further details are given in Appendix A. 2. With the numerical input for various low- 
energy constants chosen as in Ref . [8] , we find that the corrections to the isospin relations 
are very small, 

= (-0.7 ± 1.5)% , A 5 = 0.035 ± 0.003 GeV" 4 . (5.2) 

The numerical analysis shows that the corrections to a are below 1% and thus very small, 
even for isospin breaking corrections. It is interesting to note that the modification induced 
by A„ is largely counterbalanced by the modification due to A_^. Even though the modi- 
fications Eq. (5.1) only affect the rj — > 37r° tree-level couplings, these in principle also enter 
the charged channel via (inelastic) rescattering effects, however these shifts are too small 
to register. The corresponding values, which also constitute our final results, are collected 
in the final column of Table 5. 

After analyzing the isospin-breaking contributions we can now compare our final results 
for the charged Dalitz plot parameters with several other theoretical determinations and 
experimental findings in Table 6. 

We receive mixed results for the different Dalitz plot parameters. While d is in good 
agreement with experiment, a shows deviations of about 10% to the results from Layter and 
most notably from the precision measurement of the KLOE collaboration, which - due to 
the relatively small errors - exceeds even very generous confidence levels. Our result is more 
or less compatible with the 0(p 6 ) ChPT result. The dispersive calculation is somewhat 
closer to experiment, but no error range is given for us to compare with. The situation is 
even worse with b, where the deviation between our result and the KLOE measurement is 
rather alarming. The dispersive analysis indicates that even higher-order effects might be 
somewhat important in the determination of a and b, however it cannot account for the 
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Table 7: Parameterization of Dalitz plot parameters in terms of tree input parameters. 

discrepancy we find for b. A main source of uncertainty that we have not addressed so far 
is the tree- level input, which could receive rather large contributions from matching to the 
chiral amplitude at 0{p % ). It is possible that the deviation in a can be accounted for by 
such a matching prescription. There is no indication, however, that this is also the case 
for b. This issue is put under tense scrutiny in the next section. The results obtained in 
that discussion question to some extent the consistency between the charged and neutral 
Dalitz plot measurements. - - Our result for the cubic parameter, / = 0.083 ± 0.019, is 
reasonably compatible with the KLOE determination, / = 0.14 ± 0.01 ± 0.02 [27]. 

6. Relating charged and neutral Dalitz plot parameters 
6.1 Isospin limit Q n = Q c 

Up to this point we have only discussed uncertainties due to the effective range parameteri- 
zation in the final-state interactions. A by far greater source of uncertainty is the tree- level 
input to our calculation, i.e. the matching to the ChPT one-loop amplitude, which we 
deem responsible for most of the remaining deviation from the experimental results. At 
higher orders (chiral C(p 6 )), these tree parameters will receive chiral SU(3) corrections, or 
renormalizations of 0(m s ), which certainly are potentially large. In order to document our 
findings beyond the matching to the chiral one-loop amplitude, we provide a direct param- 
eterization of the various Dalitz plot parameters in terms of these input (tree) parameters. 
For this purpose, we first revert back to the case Q n = Q c as expressions become much 
simpler in this limit. Table 7 shows the coefficients of the respective input parameters. 
The entries are to be understood as follows: e.g., the second line means that the value for 
a including final-state interactions is determined by the tree input according to 

a = 0.0202 - 0.4228i + (1.0092 - 0.1902i)a tree 

- (0.0393 + 0.0182i)6 trcc - (0.0200 + 0.0378i)c? ree . (6.1) 
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All numerical coefficients are determined by tttt scattering alone. They are again averaged 
from the four different results (two loops plus unitarized, with ACGL and KPY parameters 
used as input) as in Table 4. The error range of Dalitz plot parameters calculated with 
this parameterization may be taken from the last column of Table 4. Furthermore, we only 
show the relations linear in the tree parameters (that is, no terms of quadratic etc. order), 
which are the by far dominant contributions. 

As we will now show, Table 7 can be used to construct an explicit relation between 
charged and neutral channel Dalitz plot parameters. From Eq. (2.11) one can derive the 
following relation (again, we only consider R n / R c or Q n ^ Q c in the overall normalization 
for the moment): 

* = ^(d + b-\a\ 2 ) , (6.2) 

and consequently (cf. Ref. [6]) 




which turns into an equality only for Im(a) = 0. The obvious question arises: as Im(a) 
is generated by final-state interactions but in turn depends on the Dalitz plot parame- 
ters, can we quantify the equality in Eq. (6.3) in such a way that we obtain a testable 
consistency relation between the experimental observables a, a, b, and d, independent of 
any (potentially insufficiently accurate) ChPT input? The answer is yes - precisely by 
using the information contained in Table 7. We consider Eq. (6.1) and first note that, to 
very good accuracy, the contributions from 6 tree and d tree can be neglected: with these 
parameters matched as previously, we have a trcc ps —0.656, 6 tree « —0.017, d tree ps 0.037, 
which is sufficient to demonstrate that 6 tree and d tiee are suppressed compared to a tree by 
at least one order of magnitude, irrespective of potential higher-order corrections. (The 
neglected terms are retained explicitly in the following Sect. 6.2, compare Eq. (6.9), which 
fully justifies their omission.) So via a tree in Eq. (6.1), Im(a) can be solved for Re(a) = a/2, 
and we find 

a = W c v + d ~ t) " + C2 ° )2 ' Cl = °' 050 ± °' 005 ' C2 = 0,225 ± °' 003 ' (6 ' 4) 

We wish to emphasize once more that the values for depend solely on tttt rescattering 
effects and are independent on any chiral one-loop input. The most precise determinations 
of the charged Dalitz plot parameters come from the KLOE experiment [27], see Table 6. 
Inserting their numbers for a, b, and d into Eq. (6.4), we find 

a K LOE,NREFT = -0.062 ± 0.003(stat)j£g{£(syst) ± 0.003(^tt) , (6.5) 

where the statistical and systematic errors are calculated from the respective uncertainties 
and their correlations in Ref. [27], and the last error is the uncertainty inherent in our 
assessment of final-state interactions in Eq. (6.4). This result disagrees rather strongly with 
the world average of a = —0.0317 ± 0.0016 [42] as well as KLOE's own direct experimental 
finding a = -0.0301 ± 0.0035±g;gg§§ [12]. 
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This observation seems to be at odds with a result presented in Ref. [27], where a sep- 
arate fit has been performed using an alternative parameterization [51], which incorporates 
final-state tttt rescattering based on a strict AI = 1 rule and allows to extract a therefrom. 
The result of that alternative fit is 

aA/=i,ex P = -0.038 ± 0.003(stat)l°;°J|(syst) , (6.6) 

and thus seems to be in very reasonable agreement with the direct determination of a. 
However, the parameterization from Ref. [51] is based on chiral one-loop phases or imag- 
inary parts, hence leading-order rescattering with 0(p 2 ) tttt vertices. If we reduce our 
rescattering formalism to that order (and also set Q n = Q c ), we find for the coefficients in 
Eq. (6.4) Ci = 0.021, C2 = 0.188 instead, and as a result 

a Ai= 1 ,nreft = -0.042 ± 0.002(stat)l°;^(syst) , (6.7) 

in satisfactory agreement with Eq. (6.6) within errors (which stem from the Dalitz plot 
input exclusively). We therefore understand why the rescattering formalism employed in 
Ref. [51] leads to a seemingly consistent result for a; however, the large impact of higher 
orders in the effective range parameters renders this procedure unreliable. Employing a 
more precise parameterization for tttt final-state interactions, responsible for the imaginary 
parts necessary for the relations Eqs. (6.3) and (6.4), shows that there seems to be a 
significant tension between the available experimental results for charged and neutral Dalitz 
plot parameters. 

6.2 Isospin-breaking corrections due to Q n 7^ Q c 

We now study isospin-breaking corrections to the above relations due to kinematic effects 
stemming from Q n 7^ Q c . Following the results of Sect. 5, all other effects are certainly 
included in the uncertainties. If we denote the charged Dalitz plot parameters as calculated 
in Sect. 4 by Oi so , h so , and so forth, the "real" ones a, b, . . . as deduced from Eq. (2.9) are 
related to the former according to 

a iso = a + 5 (2b - a 2 ) +0{5 2 ) , 
b iso = b + S(3f -ab) + 0(5 2 ) , 

d iso = d + 5(g-ad) + 0(5 2 ) , (6.8) 

where 5 = Q n /Qc — 1 ~ 0.069, and only / and g do not receive corrections as long as we 
disregard Dalitz plot parameters of 0(e 8 ). The corrections oc 5 produce large shifts (as dis- 
cussed in Sect. 5), so that one may wonder whether the relation Eq. (6.4) between charged 
and neutral Dalitz plot parameters may also receive large corrections. To investigate this, 
we have to amend Eq. (6.4) in two respects: 

1. incorporate the isospin-breaking shifts due to Eq. (6.8) in the terms b + d — a 2 /4; 

2. improve the parameterization of Im(a) to include C(e 4 ) effects (proportional to 6, 
d, and a 2 neglected before) in order to consistently incorporate the shifts due to 
Eq. (6.8) in the contribution to a stemming from the imaginary part of a. 
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Although significantly more complicated in result, the manipulations are much the same 
as before, relying on Table 7. The improved result is of the form 



O 2 f a 2 



a =wA b+d - 4 



- Ci [l + C2G + (Ca - ^2) a 2 + (C 4 + 2<5C 2 ) b + ( 5 d\ , Ci = 0.050 ± 0.005 , 

C 2 = 0.223 ± 0.003 , Cs = -0.008 ± 0.001 , ( 4 = 0.030 ± 0.004 , ( 5 = 0.051 ± 0.001 . 

(6.9) 

It turns out that the more refined description of Im(a) in Eq. (6.9) and therefore the 
complicated piece in the relation between charged and neutral Dalitz plot parameters 
changes the result only minimally: it shifts a by a mere +0.001. The term oc 5 in the 
first line of Eq. (6.9) is a bit more difficult to evaluate, as it involves large cancellations 
between the various contributions. This becomes evident, when analyzing the dependence 
of that term on the specific value for /. Varying / from the experimental result to our 
determination alone shifts the contribution of the 8-term from +0.002 to —0.002. 

For our final result we resort to the KLOE parameters again and use correlated errors, 
except for g, where there is neither a measurement nor a determination of its correlation 
coefficients to be found in the literature. We simply choose to vary it independently 
between zero and the result of our calculation in Table 5. However, despite these generous 
variations, the total effect of these additional contributions proportional to S is still so 
small that it hardly shows in the overall uncertainty. Our final result is 

a K LOE,NREFT = -0.059 ± 0.003(stat)t°;^(syst) ± 0.003(^tt) . (6.10) 

The overall correction to Eq. (6.5) turns out to be small and we are still left with a 
significant disagreement between charged and neutral channel. Comparing the charged 
Dalitz plot parameters entering Eq. (6.4), we see that the main disagreement is due to the 
parameter b, which is strongly over-predicted in our analysis: we find b = 0.308 ± 0.023 
to be compared with 6 K LOE = 0.124 ± 0.006 ± 0.010. (Of course, the NREFT results are 
consistent within themselves: inserting our values for a, b, d into the relation Eq. (6.4) 
reproduces our result for a.) We also mention that there is some non- negligible variation 
between the KLOE results for the charged Dalitz plot parameters and several older, less 
precise measurements [49,50,52]; a re-measurement of these quantities by some of the 
modern high-precision experiments would therefore be very welcome. 

The relation between a and the charged Dalitz plot parameters is further illustrated in 
Fig. 4. Due to the smallness of the higher-order corrections in Eq. (6.9), it suffices to use the 
simplified representation Eq. (6.4). As a = 0.0317 ± 0.0016 is experimentally agreed upon 
to very high precision, and as our result for d agrees well with the KLOE determination 
d = 0.057 ± 0.006lio oi6' we ma y * a ^ e these two experimental results for granted, such 
that Eq. (6.4) provides a relation between a and b. This constraint in the a — b plane is 
shown in Fig. 4. The solid grey area shows the allowed range for 6 as a function of a 
according to Eq. (6.4), whereas the hatched grey area shows the same relation for Q\ = 0, 
i.e. fully neglecting the imaginary part of the amplitude, or Im(a) = 0. While the NREFT 
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prediction for a and b falls nicely into the allowed band (the agreement here looks even 
better than in the direct comparison to a as the band also reflects the experimental error 
in d), the KLOE determination of both is consistent with a vanishing imaginary part. 
In our framework these latter values cannot be brought into agreement with a consistent 
implementation of final-state interactions. 

Comparing our calculation to the dispersion-theoretical analysis of Ref. [22], there are 
indications that the discrepancy we find may be slightly over-predicted: in the terminology 
of the iterative solution determined there, our two-loop calculation cannot be expected to 
be better than the second iteration of the dispersive amplitude. Ref. [22] shows that while 
the real part of the amplitude has converged to the final result almost perfectly, there 
are still non-negligible corrections in the imaginary part beyond that, i.e. in terms of our 
representation at (irreducible) three loops and higher. Whether those corrections in the 
imaginary part that precisely constitute the additional terms in the relations Eqs. (6.4), 
(6.9) are sufficient to reduce the discrepancy between charged and neutral Dalitz plot 
parameter measurements remains to be seen. 



7. Partial widths and the ratio r 

In this work, we have concentrated almost exclusively on the energy dependence of the 
two 7] — > 3ir Dalitz plot distributions, mainly as encoded in the Dalitz plot parameters. It 
is rather obvious in particular from Table 3 that the overall normalization of the ampli- 
tudes is not improved in our formalism compared to what we match our parameters to, 
here ChPT at 0(p 4 ); indeed, the overall rates are even slightly smaller. To be concrete, 
integrating the Dalitz plot distributions of our amplitudes including all isospin-breaking 
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effects (corresponding to the last column of Table 5), we find 

T(r] -»• 3^°) = [201 ± 3(^tt) ± 6(A_#-)] eV , 
r(r? -> 7r + 7r _ 7r ) = [144 ± 2(vr7r)] eV . (7.1) 

Several remarks are in order here. First, as we have pointed out earlier, for our normal- 
ization we use Q = 24.2 as given by means of Dashen's theorem, which leads to a very 
small width. Changing the value to Q = 22.3 [22], say, immediately increases the widths 
by nearly 40%. Second, Ref. [6] finds that next-to-next-to-leading order chiral corrections 
increase the width by nearly 70%, thus bringing it a lot closer to the experimental value 
of about T(r] — > 7T + 7r~7r°) ps (296 ± 16) eV [42]. We wish to emphasize once more that this 
failure to reproduce the chiral enhancements in the width in the non-relativistic framework 
does not invalidate our predictions for the Dalitz plot parameters: the power counting ar- 
gument of Sect. 3.4 explains why we catch the important rescattering effects in particular 
for the higher-order energy dependence, but not in the overall normalization. The 77 — > 37r 
tree-level coupling constants that receive sizeable quark-mass renormalization effects nicely 
factor out of the complete (tree plus loop) amplitudes and play no role in the calculation 
of the Dalitz plot parameters. As a third remark, the errors shown in Eq. (7.1) do not 
at all reflect these uncertainties from our matching procedure, but purely the one due 
to -K7T final-state interactions (determined as in the previous sections), and in the case of 
T(rj — > 3tt°) due to the uncertainty in A^-, see Eq. (5.1). 

Despite all the above-mentioned deficits in a calculation of the decay widths, the ratio 
of neutral-to-charged partial widths r should be predicted much more reliably, as the 
normalization of the amplitude (in the isospin limit) drops out. In particular, here we may 
expect a somewhat heightened importance of isospin-breaking corrections [8]. We find 

r = ->• 37r °) = 1 40 ± o.01(7T7r) ± 0.04(A^) , (7.2) 

T(r] -»• vr+vr-vr ) v ' y ' 

in agreement with the experimental finding r = 1.43 ± 0.02 [42]. We note that the depen- 
dence on 7T7T rescattering in Eq. (7.2) is very small, our error is dominated by the 1.5% 
uncertainty in A^. Equation (7.2) is extremely accurately reproduced by just integrating 
the phenomenological Dalitz plot distribution, with our values for the Dalitz plot param- 
eters from Table 5 (last column) instead of the exact amplitudes: obviously r is affected 
by cusps in the neutral channel or yet-higher-order Dalitz plot parameters at or below the 
permille level. We can therefore easily derive the dependence of r on the parameters a, 
b, d, . . . , making use of the relation of a to these in Eq. (6.3) and neglecting pieces that 
affect r at the permille level (e.g. the terms oc j3, 7 in the neutral rate), and find 

r = 1.485 (l- 0.029 a- 0.061 a 2 +0.024 6+0.032^+0.008/ -0.014 ( 1 + ^£ ) . (7.3) 



Errors on this result are to be taken from Eq. (7.2). The various numerical coefficients are 
given by tttt phase shifts and phase space integration only. This demonstrates to very good 
approximation that r does not depend on the normalization and thus possibly sizeable 
quark mass renormalization effects. 
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8. Summary and conclusion 

In this article we have analyzed rescattering effects in rj — > 3-ir decays by means of the modi- 
fied non-relativistic effective field theory framework. The main findings of our investigation 
can be summarized as follows: 

1. NREFT provides a simple and transparent representation of the amplitude to two 
loops, including higher-order isospin breaking. In order to estimate higher-order 
loop effects we have furthermore applied a simplified unitarization prescription. The 
amplitude thus obtained is - at the very least - fully competitive with the chiral 
expansion at next-to-next-to-leading order. The coupling constants involved have 
been matched to phenomenological tttt scattering threshold parameters and, in the 
case of the r\ — > 3ir tree-level couplings, to ChPT at 0(p A ). 

2. One- and two-loop contributions to the Dalitz plot parameters are in general of 
the same size, an observation which is predicted by non-relativistic power counting 
arguments. Irreducible two-loop graphs are generally suppressed, while derivative 
couplings at two-loop level are essential to find the correct sign for the rj — > 3tt° slope 
parameter a. Higher-order effects beyond two loops were shown to be relatively small, 
but not negligible. 

3. While our results for the Dalitz plot parameters are in qualitative agreement with 
previous dispersive results, we can provide an explanation for the apparent failure of 
two-loop ChPT to reproduce a: the treatment of tttt final-state interactions is still not 
sufficiently accurate at that order. We can identify one specific diagram, the double 
rescattering graph with tttt vertices beyond leading order, as being responsible for at 
least half of the discrepancy between the 0{p & ) prediction for a and the experimental 
value. These effects are of chiral order p s and higher, but included in the NREFT 
two-loop representation. 

4. Apart from normalization effects and subtleties in the definition of the center of the 
Dalitz plot in the charged decay channel, higher-order isospin-breaking corrections 
on the Dalitz plot parameters are very small. 

5. Our final result for neutral Dalitz slope parameter, 

a = -0.025 ±0.005 , (8.1) 

is compared in Fig. 5 to several other determinations. It is considerably closer to 
the experimental world average a = —0.0317 ± 0.0016 [42] than previous theoretical 
approaches. Notice though that Eq. (8.1) does not take uncertainties stemming from 
matching to ChPT at 0(p 4 ) into account, which we expect to be non-negligible. 

6. Our results for the charged Dalitz plot parameters show somewhat larger deviations 
from the currently most accurate measurement by the KLOE collaboration. By 
relating charged and neutral decay channel via the AI = 1 rule we find indications 
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Figure 5: Comparison of values for the slope parameter a. Top: theoretical predictions. Bottom: 
experimental determinations. The grey shaded area is the particle data group average [42]. 

for a significant tension between the Dalitz plot parameters of both channels, which 
is solely due to final-state interactions. A re-measurement of the charged Dalitz 
plot parameters by high-precision experiments [53,54], or even preferably access to 
improved full Dalitz plot distributions, is thus highly desirable. 

7. While the partial widths calculated in our framework do not improve upon the chiral 
one-loop prediction we match to (due to the absence of further quark-mass renormal- 
ization effects not captured in our framework), we can give a value for the ratio of 
neutral-to-charged partial widths unaffected by this deficit, r = 1.40 ± 0.04, where 
the error is dominated by isospin-breaking effects. 

Possible future improvements on the theoretical approach include matching to C(p 6 ) 
ChPT in order to constrain the tree-level Dalitz plot couplings more tightly. Furthermore, 
it will be extremely useful to match the non-relativistic representation to the upcoming 
dispersive analysis [22] in order to obtain a reliable description of the whole physical Dalitz 
plot: in this way one can include elastic tttt rescattering to all orders, and at the same time 
implement in particular non-analytic effects (cusps) at or near the boundaries of the Dalitz 
plot due to isospin-breaking up to next-to-next-to- leading order. This combination should 
then also provide the best-possible representation of the decay amplitude for a precision 
extraction of the quark mass ratio Q. 
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A. Isospin-breaking corrections 

A.l Isospin-breaking corrections to 7T7t scattering 

To calculate isospin-breaking corrections to the matching relations Eq. (3.7), we expand 
the ChPT amplitudes for all channels with electromagnetic corrections included [47, 48] 
around threshold. These contain virtual-photon exchange and real-photon radiation in the 
form of bremsstrahlung. For the definition of a reasonable (regular) threshold expansion, 
at first the divergent Coulomb pole contribution has to be subtracted. As the Coulomb 
pole emerges equally in the vertex correction diagram in both the NREFT and the ChPT 
calculation due to the same infrared properties of both theories, in a matching between 
them, this part drops out anyway. The determination of the scattering lengths is then 
relatively straightforward and has already been performed in the above references. Note, 
however, that in contrast to Refs. [47, 48] we expand around an isospin limit defined in 
terms of the charged pion mass. Thus, non-analytic terms oc y/ A n arise in the expansion of 
the 7r°7T — > 7r°7r° channel, which are due to a cusp structure at the charged pion threshold 

and cancel the corresponding contribution in the expansion of (s) in Eq. (3.3), once the 

correct matching is performed. We display the corrections in the form Cj = Q+ACj, where 
Ci denotes the corresponding coupling in the isospin limit (the tt + tt + — > it + it + channel is 
not needed in the present analysis, we just give it for completeness). For the combinations 
of S-wave scattering lengths, we find 

M 2 r A p 2 A 
AC » = if I" 5? + 32^ C + 32^%< 13 - 16h - + 6h - *>} ' 

M 2 c A p 2 A i 

AC - = f { M - 1^ (2 ° - + l^J (3 + 2h + 4,;) } ■ {AA) 

We note that ACoo is indeed free of non-analytic terms in A^: the analytic structure of 
ChPT and the non-relativistic representation near threshold is the same, as it must. 

The definition of effective ranges and P-wave scattering lengths is not a priori clear, 
since one has to deal with infrared divergences in the ChPT amplitudes. In calculations 
of, say, cross sections these divergences, which arise from virtual-photon corrections, can- 
cel with corresponding divergences from real-photon radiation (bremsstrahlung). How- 
ever, when matching the non-relativistic framework to ChPT, the explicit inclusion of 
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bremsstrahlung is not necessary, since the virtual-photon diagrams exhibit the same in- 
frared behavior and thus contain the same divergences (see Ref. [31]). On a rather 
technical note, the infrared divergences were calculated in dimensional regularization in 
Ref. [31], while the ChPT calculations [47,48] use a finite photon mass m 7 as infrared 
regulator. The transition between both regularization schemes can be made by replacing 
log(m 2 /M 2 ) — 32-7t 2 Air — 1- The infrared divergences then cancel, rendering the match- 
ing relations finite. We wish to emphasize that the physical reason for this cancellation is 
again the identical infrared behavior of both theories. 

For the S-wave effective ranges, defining in analogy with the above Dj = Di + AD. 
we find the following corrections to the matching relations: 



ii 



AD ++ = j} {- 5^( 676 + 9 < c+ - - + iss^} • < A - 2 > 

while for the two P-wave scattering lengths, we have (with Ei = Ei + AE{) 



-'— r - 28 « + - K++ )-i^k}- < A - 3 » 



AE+ 



Fl 1 128vr 2 



where the following abbreviations have been used for combinations of electromagnetic SU(2) 
low-energy constants fcj and Z = A n /(2e 2 F%): 

= ( 3 + ^) ~ kl - ^k 2 - 3k 3 ~ 4^4 , 



4Z \ - 32Z - - 
/cf u = ^3 + -)*! + —h + 3h + 4Zk 4 



-±0 

-1 

ICf = 8Zk 2 + 3k 3 + AZk A - 2(1 + 8Z)k % - (1 - 8Z)k 8 , 

K+- = (3 + ^) fci - -^fc 2 - 9fc 3 + 4Zfc 4 + 4(1 + 8Z)fc 6 + 2(1 - 8Z)k 8 , 

K + - - = 2(3 + ^) h + - 18fc 3 + 24Zh ■ (A.4) 

We refrain from calculating corrections to the shape parameters, since their intrinsic, 
isospin-symmetric error is much larger than what can be expected from isospin breaking. 

For the numerical evaluation we express the low-energy constants l\ and l 2 in terms of 
7T7T D-wave scattering lengths [55] , for which we use the numerical values [39] 

al = 1.75 ± 0.03 x 10~ 3 M~ 4 , a\ = 0.170 ± 0.013 x 10~ 3 M~ 4 . (A.5) 
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channel 


ACi/d x 10~ 2 


A A/ A x 10" 2 


AEi/Ei x 10~ 2 


00 


-7.3 ±0.2 


-3.3 ±0.4 




X 


2.5 ±0.6 


0.1 ±0.6 




+0 


-5.2 ±0.8 


1.9 ±0.8 


0.4 ±0.6 


+- 


6.1 ±0.5 


-0.2 ±0.3 


0.5 ±0.4 



Table 8: Corrections to the matching relations relative to the phenomenological values. 



This way a® and a 2 , can be independently varied according to their uncertainty, whereas l\ 
and I2 are correlated. For l 3 we propose 1% = 3.1 ± 0.5 as a sensible mean value from lattice 
simulations (see Ref. [56] for individual results of the various groups). The constant U is 
extracted from the scalar radius of the pion [39], I4 = 4.4 ± 0.2 . 

For the electromagnetic SU(2) low-energy constants k\ we use the values given in 
Ref. [57]. The authors of this work have matched the two-flavor low-energy constants to 
their SU(3) counterparts, using numerical estimates from Refs. [58,59]. We convert the 
values k\ given at the mass of the p, M p = 0.77 GeV, in Ref. [57] to scale-independent 
constants according to the standard prescription, 

327r 2 M 2 
k = — fc[(M p )-log-|, (A.6) 

"i lvl p 

where the <Tj are the corresponding f3- functions to be found in Ref. [47]. Numerically this 
results in 

h = 1.66 , k 2 = 4.08 , k 3 = 2.28 , fc 4 = 3.69 , k 6 = 4.08 , k$ = 4.06 . (A.7) 

The uncertainties on the k\ are estimated analogously to Ref. [8] by their logarithmic scale 
variation, 

K -> K ± ^ , (A.8) 

which for the fcj translates to k\ — >■ ki ± 2 . The errors on the quantities /C 00 , JCf°, K+~ , /C++ 
are then calculated in a correlated fashion (i.e. ±2 or —2 for all ki). 

The numerical corrections are displayed in Table 8. We find that the corrections at 
one-loop order are very small. The main contributions to the Q stem from the tree-level 
correction factor. 

A. 2 Corrections to the AI = 1 rule 

At leading order p 2 in ChPT and up to next-to- leading order in the isospin-breaking param- 
eters m u — rrid and e 2 , the amplitudes for the charged and the neutral decay were already 
quoted in Eq. (2.6). We also hinted at the fact that in order to define the deviations from 
the AI = 1 relation, Eq. (5.1), it is useful to expand the decay amplitudes for both channels 
around the point S3 = s n , s± = s 2 as shown in Eq. (2.8), so that is going to be of chiral 
order p 4 . 
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The decay amplitudes at 0(p 4 ) in ChPT and at 0{m^ — m u , e 2 , {m^ — m u )e 2 ) in isospin 
breaking are given explicitly in Ref. [8]. With minimal modifications they can be shown 
to be also valid up-to-and-including 0{{m,d — rn u ) 2 ), i.e. only numerically tiny terms of 
0(e 4 ) are potentially neglected at second order in isospin breaking. In order to match the 
expanded ChPT amplitude of Ref. [8] to the polynomial part of the NREFT representation, 
the following steps have to be taken into account: 

1. The (non- analytic) imaginary parts due to pion loops in chiral and NREFT amplitude 
are identical and drop out in the matching relation. 

2. For the radiative corrections due to real and virtual photons, we have to match the 
result of Ref. [8] to an analogous NREFT representation as in Ref. [31]. As a result, 
the Coulomb pole and phase have to be subtracted from the chiral representation, 
as well as the bremsstrahlung contributions. As in the case of radiative corrections 
to 7T7T scattering infrared divergences were regulated by introducing a finite photon 
mass m 7 in Ref. [8] and have to be treated as described in Appendix A.l. 

As Ajy- is of 0(p 4 ), it is convenient to factor out the neutral normalization at leading 
order and quote the result as the ratio Ajj-/J\f n below. We find 

^ - 2e 2 { 1 —^G(s n ) + -l)- \±^(2Kl - KJ) 



M n 1 3p K nj ' 2 wny nj ' 32vr 2 V V 2 J 1-p 

+ 3(1-/9) 1-/9 ^10+^ll) 

, A„ / 29-lllp-9p 2 + 27p 3 = 

+ 3(l-p)F|\ W^P) KK{ n) " 3 

3/9(1 + 22p + 9p 2 ) - 7 + 21/9 - 495/9 2 + 243/9 3 - 

+ (1 _ 9/9 )(i + 3/9 ) Vl^J (1 _ Qp)(1 + 3p) J ™{Sn) 



8(3-/9) F 2 , 1 



1-p M2 * 16vr 2 



p z 1 — p 4p 



3(3 - 26p - p 2 ) , 53 - 357p + 351p 2 + 81p 3 

logp + 



(l-9p)(l-p) ^ 4(1 -9p) 

Q e 2 f 

I -Up- 138p 2 + 234p 3 - 1107p 4 - 

■Jmr\ s n) 



}+0(eV) , (A.9) 



+ 



1 - 9p 

7 - 102p - 504p 2 + 1926p 3 - 3375p 4 
32vr 2 (l - 9p) 



3A^ f3(3 + p)(957- 5240p- 1398p 2 - 288p 3 + 81p 4 ) - 

+ (l-p)(l + 3p) 2 F2M4\ 4096(1 + 3p) ^x(*„) 

p 2 (221 - 3612p + 32022p 2 - 32076p 3 - 2187p 4 ) 7 
+ 8 (l-9p) 3 (l + 3p) vASn) 

3p(3 - 124p + 1794p 2 - 7596p 3 + 9315p 4 ) - 

(l-9p) 3 (l + 3p) ^ (Sn) 
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+ 32768^(1 -9 P )4 243 - 39737 " + 540m » 2 - 729333 " 3 
+ 3630825/? 4 - 1810107p 5 + 85293p 6 + 59049//) 

p 2 (37 + 237» - 2025p 2 + 3159p 3 ) ] m . , 2N . , . 

- 128^(1 -9,)>(1-,) - '°^j+0(-=V 2 ) , (A.10) 

neglecting even higher-order terms in the isospin-breaking parameters e 2 and rtid — m u . 
Here, we have used Ap = Fk/F w — 1 (cf. Ref. [60]) and the loop functions 

G( S ) = ^#(Li(i^)-Li(i±^)+logi±^l, Li(*)= f 
Jniris) = |l - Y^ rr^ } ' Jkk{s) = {1 - cr^arccotcr^} , 

J v*( 8 ) = ^^{ 2 + l °sp[ — „ i — : ) - t 1o s 



32vr 2 1 tor V s 1-P/ s s-M2-M T 2 + y 



is the real part of the triangle loop function for the photon exchange between two 
charged pions (rescaled by a factor of M 2 ) with the Coulomb pole subtracted, involving 
Spence's function Li(z), and J a b(s) are the usual finite and scale-independent parts of the 
corresponding two-meson loop functions. For the definition of the (renormalized) strong 
and electromagnetic SU(3) low-energy constants Lg and in terms of chiral Lagrangians 
(not to be confused with the tree- level couplings Lj, Ki of the non-relativistic theory), see 
Refs. [60,61]. We have made extensive use of the Gell-Mann-Okubo relation to simplify 
the results Eqs. (A. 9) and (A. 10). Furthermore, the kinematic expansion to second order 
plus the expansion in isospin-breaking parameters leads to derivatives up to third order of 
these loop functions, which have been rewritten in terms of the loop functions themselves. 

Both results Eqs. (A. 9) and (A. 10) are divergence-free and independent of the scale [i. 
The scale-independence of A^ is explicitly seen, that of A^ can be found by using the scale 
variation of the electromagnetic constants K\ as given in Ref. [61]. Both corrections turn 
out to be completely of electromagnetic origin. For the numerical evaluation we make use of 
the same estimates and variations of the low-energy constants as described in Appendix A.l 
and explained in more detail in Ref. [8]; their uncertainties completely dominate the error 
on AjQ-/M n . As A„ is free of low-energy constants at this order, it is a pure loop effect and 
a prediction in terms of well-known parameters. Here we quote an uncertainty solely due 
to the use of the Gell-Mann-Okubo relation for the masses, using either the mass of the rj 
directly, or the same expressed in terms of pion and kaon masses. We consider the error 
thus obtained rather underestimated. In total, we find 

^f- = (-0.7 ± 1.5)% , Aa = 0.035 ± 0.003 GeV -4 . (A.12) 
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B. NREFT representation including isospin breaking 
B.l rj — > 3tv amplitudes up to two loops 

For the representation of the 77 — > 3ir decay amplitudes at one-loop order, we find (see also 
Ref. [30]) 

M 1 n - loop (s 1 ,s 2 ,s 3 ) = [Coo(si)K(s 1 )J 00 (s 1 ) + 2C x (s 1 )L(s 1 )J^(s 1 ) + ( Sl o s 2 ) 
+ (*i s 3 )} , 

A^c" loop (si,S2,s 3 ) = C X (S 3 )K(8 3 )JOO{S3) + 2C + -{8 3 )L(S 3 )J+-{S 3 ) 

+ { [2C+o(si)L'( Sl ) - E+o( Sl ,S2, ss)^(si)] J+o(«i) + (*i ^ s 2 )} • (B.l) 

The topologies at two-loop order are shown in Fig. 2. For the two-loop amplitudes at the 
order discussed in Sect. 3.4, we obtain 

M 2-loops = S2j S3 ) + M B( Su S2j S3) + (si S2 ) + (si o S3 )} , 

_ M 2-loop S = _ M A (si)S2)S3)+ _ M B (si)S2)S3) ; (R2) 

where 

A< = 2K Coo(^)C'oo(si)Fo(M 7r o,M 7r o,M 7r o,M 7r o, Sl ) 

- 4M v K D 00 C 00 ( Sl )^F^\M n o,M n o,M n0 ,M n0 , Sl ) 



+ 8 



L'|' ( Sl )C +0 (5+") 
A n /si + 2Qf 



V 2Q" 



p? ) L £ +0 (ai)F (M^o , M w , M n , M n , Sl ) 



rl 



-LiC+o^-) + 2M r) L';' ( Sl )£»+o 



4M 7r o ' 



C-( s i)S F o - M ^ M ^ M ^ si) 



(OS] 



+ 8M r ,L 1 £» +0 C' a; ( Sl ) 7 ^ J P (2) (M 7r o, M,,, M,,-, M n , Sl ) 



+ 4L"( Sl )C a; (s? )C o(s 1 )Fo(M 7r , M ff , M n o,M n o, Sl ) 

Q? *<i>, 



+ 4 



LiC x {S?) -2M v L"( Sl )D x Coo(si)^F^ , (M 7n M n ,M n o,M n o,s 1 ) 



Qi 



(B.3) 



At* = K(si)C 00 (si) 2 J 2 ( Sl ) + 2[L( S i)C x (si)Coo(5i) + K(si)C 2 ( Sl )] J 00 (si) J+_(si) 

+ 4L( Sl )C + _(s 1 )C a; ( Sl )J 2 _( Sl ) , (B.4) 



At 



L , ;' ( S +)C + o(S+ )(c +0 (si) - ^0(^1^2,53)) 
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A, / Sl + 2Qf-A, 



4M 7r o 



2Q\ 



- p?) L E +0 C +0 (si)] F + (M w , M^o , , M^o , «i ) 
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+ 4^ LojB+ ° C+o(Sl) ^f 



+ 4 
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p. 4 
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+ 8A^LiD+o£; + o(si, 8 2 , s 3 ) 7 ^F^>{M 7T ,M 7T o,M n , M n o, Sl ) 

Wi) 
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Qi 1 



- 2M v L"(s^)D + ) E +0 ( Sl , 82, s 3 )^] F^(M„, M„, M n o,M n , 8l ) 

+ 16M r) L 1 D + -E + o{s 1 , 82 , 33)7^2 F+\M V , M n , M n o,M n , 8l ) 
+ 2KqC x (s°+) (c +0 (si) + E- Q ( Sl , 8 2 , s 3 )^F + (M n o, M n o, M w o, M n , si) 

Q? 

Q? 

p(i) 



4if n 



M V D X (c +0 (si) + £+ (si, s 2 , s 3 )) ^ + C :r (s° + ) J E; + o(si, s 2 , a 3 ) 
x (M w o , M w o ,M„o,M„,8i) 
+ 8M v K D x E +0 (s 1 ,82, s 3 )^F^(M n o, M n o, M n o, M W) si) + ( Sl o s 2 )| 
+ 2K C , 00 ( S T)C a; ( S 3)Fo(M 7r o, M^o, M u o, M,o, a 3 ) 

q2 

- 47^^0^00^(53)^ F (1) (M w o, M w o, M w o, M,o, s 3 ) 
V 3 

+ ±L"(s 3 )C x (^)C x {s 3 )F Q {M n , M n , M n o,M n o,s 3 ) 

LiC x (%>) - 2M V L"( S3 )D X ] C x (s 3 )^F^(M n , M n , M^,M^,s 3 ) 

Q 3 



+ 4 



Qi 



8M v L 1 D x C x ( 83 )-^FP(M n ,M n , M n o, M w o, s 3 ) 
W 3 ) 
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+ 8 



L'^(s 3 )C +0 (~st-) 
A n /s 3 + 2Qi 



4M 
1 



7,( 



2Q% 



pf) L E +0 C+_ (53)^0 (M w , M n0 , M w , M w , s 3 ) 



LiC +0 (4-) + 2M V L'% (s 3 )D +Q 



12 



A 



7—LnE. 



C+- (as) ^|^ (1) , M,,., M,,., s 3 ) 

+ 8M V L 1 D +0 C+- (s 3 ) F (2) ( M - , M w , M n , s 3 ) , 



(B.5) 



Mf = {4L'( Sl )C2 ( Sl )4 ( Sl ) + ( Sl o s 2 )} + K( S3 )C 00 (s 3 )C x (s 3 )Ji (s 3 ) 

+ 2[l(s 3 )C 2 x (s 3 ) + ^(53)^(53)^(53)] J+-(s 3 ) Joo(s 3 ) + 4L(s 3 )C 2 + _(s 3 )Jl_(s 3 ) . 

(B.6) 



We have used the following abbreviations: 



Jab{Si) 



\( Sl ,MlM*) 



iQab{Sjj 

„thr\ 1 77 /„ „thr N \ 



C n (si) — C n + D n (^Si — s\ r ) + F n (s{ — r ) , 

5~f = M c 2 + M? - 5, + ^ (5, + 2Q? - M 2 + M 2 ) , 



^+0(S1,S2,S3) = -^+0 
£ + _(si,S2,S 3 ) = E + _ 



^W( S1(S3 _ S2) + M M|-M ?) ) 
(si(±A^))(s 3 -s 2 + A 7r ) 



2M V Q\ 
sz{si - s 2 ) 
2M v Ql ' 



-A, 



K( Si ) = ^o + ^i (p?-M^ + 2 ^-M,o + ^ 1 



Q? 



Qi 



0? 



L( Si ) = L + Lx (p° - M w o) + L 2 (p° - M^o) + L 3 ^- (l - — 



4M 2 

Si 

4M 2 \ 



L(si) = Li + 2L 2 



+ 2L 3 



0? 



+ £3 



s 



K'M^) =L + Ll (\(M v -M a - M b ) - s -^l) 



(B.7) 
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In the notation, it is understood that the shape parameter term F n is omitted in the 
polynomials C n (sf d ) inside the "genuine" two- loop graphs. There is a subtlety with regard 
to the neutral Dalitz plot couplings in the irreducible two-loop graphs: since we included 
these couplings only up to 0(e 2 ), the AI = 1 rule is only fulfilled up that same order. 
Eq. (3.2), however, is valid up to C(e 4 ), so that we have to replace Kq — > K$ = — (3Lq + 
L\Q n ) in M-n and above or simply Kq = J\f n (l — ^bM^Q^). The numerical effects of 

this replacement are small. 

(k) 

Fi(...;s), (. ..;s), k = 1, 2, 3, stand for the integral representations F(...;s), 
F^(. . . ; s), evaluated at Q? = A(M 2 , M 2 , Si)/AM 2 , with i = 1, 2, 3. The analytic expres- 
sion for these two-loop functions read 



30 2 

F(M a , M b , M c , M d , s) =J\f 2Af 1 + Bf - -£-(B h + 2C f ) 



Ws 



+ K(X 3 f 3 + X 2 f 2 + X ± h + X fo) 
N 



F^(M a ,M b ,M c ,M d ,s) = J -^(l + 5)[(10A-B)f 1 + (5B-2C)fo}+0(e 4 ) , 
F^(M a ,M b ,M c ,M d ,s) = ^ 



^ (2A 2 f 3 + 3ABf 2 + (B 2 - 2AC)h + BCf ) 



+ 



(1 + 5)' 



(a/ 3 + (B- 2A)f 2 + (4A -2B + C)h + 2(B - C)/ ) 



+ 0(e 4 ) , 



F^(M a ,M b ,M c ,M d ,s 



M(l + 8) 



16 



3 



(A 2 /4 + 2A(B - 4A)/ 3 

+ (2AC + B 2 - \2AB)h + (2BC - 8AC - 4B 2 )/i + (C 2 - 4BC)/ ) 
(1 + 5f (a/ 4 + (B - 3A)/ 3 + (3A - 35 + C)/ 2 



+ (3B - 4A - 3C)/i + (3C - 25) / 



(B.8) 



with 

/2 = 
/4 = 



■— (3B/i + C/o), / 3 



7A 2 



3(AC - 5 2 )/i - BCfo 



9A 2 



3(A5C7 + 5BC - 5B 3 )h - 5(B 2 C - C 2 )f Q 



X = BBC -RC, Xi= H(2AC + B 2 ) - R(2B - C) , X 2 = 3HAB - R(3A ~^B), 



X 3 = 2HA 
and 

N = 



2AR, H 



2 V 3a 



R 



( i + if 



(B.9) 



X 1 / 2 (s ,M 2 ,M 2 ) 



/C 



1 



+ 



2(M 2 + M 2 ) - (M a + M b ) 2 - s So -(M a -M b ) 



so 



M, 



«0 



M 2 - M 2 
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/o 

fl 

h 

Vi 
2/1,2 

A 



4(v 1 + f 2 - V 2 + h) , 
4 
3 

im, 

2 V.l + Q 2 /s 



(yi(«i - 1) + y 2 (v 2 - 1) - 2/2(^2 - 1) + h) , 
l + Q 2 /^ 



Q = Q (s) , 



Mr/) 2 



\J—yi arctan • 



1,2; f 2 = V - 2/2 arctan 



-5 T V5 2 - 4^C 



2.4 



1/2 = 1/2 (s) , 



Q 2 



--j- (M c + A^ 



O 2 

£ = - A 2 + — M 2 C 



C 



s = M 2 + Ml - 2M V Udl + 
x2 _ A(M 2 ,M 2 ,(M Q + M b ) 2 



2 , Q 2 (i + ^7 



2\ 1/2 



4M 2 



B.2 Resummed amplitudes 



~ M 2 



2 _ A(g,M 2 ,M 2 ) 



(B.10) 



In order to estimate the effects of higher-order corrections we iterate the bubble diagrams 
and the external vertex of the non-trivial two-loop graph. A diagrammatic expression of 
this iteration is shown in Fig. 3. Here we show the results including isospin violation. For 
the bubble chain a coupled-channel resummation can be performed analogously to Ref. [62]. 
We obtain 



Mn(s!,S2,S 3 



2L(s 1 )C x (s 1 )J+-(s 1 ) + K( Sl ) [Coo(si) Joo(si) - 2 x (si) J+-(«i) Joo(*i) 



1 - 2C + _(si)J + _( Sl ) - Coo(si) Joo(si) + 2 x (si) J+-(si) Joo(si) 
+ (si s 2 ) + (si o s 3 ) , 
A/fU/ x 2C +0 (s 1 )J +0 (s 1 )L / (s 1 ) E' +0 (s 1 ,s 2 ,s 3 )L(s 1 )J +0 (s 1 ) 

M C {S 1 ,S 2 ,S 3 ) = — — — — — — — r h (Si O S 2 ) 



+ 



1 - 2C +0 (si)J+o(si) 1 - E + o( Sl )J + o(si) 

2L(s 3 ) C + -(s 3 )J + -(s 3 )-2 X {s 3 )J + -(s 3 )Joo{s 3 )\ + K(s 3 )C x (s 3 )J 00 {s 3 ) 



1 - 2C + _(s 3 ) J+-(s 3 ) - C o(s 3 )Joo(s 3 ) + 2 x (s 3 )J + „(s 3 )J 00 (s 3 ) 



(B.ll) 



where 



X(sj) = C + _(si)C o(sj) - C x (si) 2 

J2, 



E' +0 (s!, s 2 , s 3 ) = E +0 + G +0 ( Si - s%) 



_^( <i(j ,_ Jj) + M «._«J ) ) 



^+o(si) 



4«+o( a i) 



E +0 + G +0 ( Si - s%) 



G 



+o 



37rfei , 



(B.12) 



and b\ is the P-wave effective range. Additionally, we performed a resummation of the 
external vertex of the non-trivial two-loop diagram. This can be achieved by replacing the 
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Figure 6: Two- loop graphs with three-particle cuts in the physical region. 



outer vertex in Eqs. (B.4) and (B.6) according to 

Coo(si) - 2x{si)J + -{si) 



Coo(si) 
C x (si) 



1 - C 00 (si)J 00 (s,i) - 2C + -{si)J + -{si) + 2x(si)J + -{si)J 00 (si) 

Cx( s i) 

1 - C o(si)J o(si) - 2C + -(si)J + -(si) + 2x(si)J+-(si)J o(si) 

C+o(si) ( ± ) E' i ^(s 1 ,S2,S 3 ) 



C + o{si)^- — , Ey(si,s 2 ,s 3 ) -»• , 

1 - 2C +0 (si)J+o(si) ^ 1 - E +0 (si)J +0 (si) 

r C + -(sj) - x(si)Joo(si) 

+ ~ {l> l-Coo(« i )Joo(s i )-2C + _( Si )J + _( 5i ) + 2x(«i)J + -(« i )Joo(«i) ' l ' J 



with 



thr\ 



(si(±A^))(g 3 -g 2 + A^) 



C. Comment on imaginary parts of two-loop diagrams 



(B.14) 



In our analysis of the non-relativistic rj — > 3tt decay amplitude, we have neglected the 
imaginary parts of the non-trivial two- loop graphs, see Fig. 6 (left). The loop function 
F(s) (in the simplified notation introduced for the equal-mass case in Sect. 4) given in 
Appendix B.l strictly speaking only corresponds to the real part of this diagram. At 
leading order in the e expansion, its imaginary part is given by 

V ' V3(32vr) 2 Ml K ' v ' 

We therefore confirm that ImF(s) = C(e 4 ), while the real part of the same diagram already 
starts at 0(e 2 ). The imaginary part is due to the three-pion cut and only arises because the 
rj is unstable, M„ > 3M n . It stems from a part of the non-relativistic loop integral in which 
one of the propagators is non-singular, and therefore yields a result very similar to that 
of the sunset graph Fig. 6 (right), which in the non-relativistic framework can only arise 
when introducing (very small) six-pion vertices [29,63]. The three-pion cut causes a non- 
vanishing imaginary part of the isospin amplitudes Mi(s) already at threshold s = 4M 2 . 

It is obvious that ImF(s) can only contribute to the amplitude's squared modulus 
at 0(a^ r7T e 5 ) via interference with one-loop terms, and is therefore naturally suppressed 
compared to the real part at two-loop order. What is less clear is its relative importance 
compared to the imaginary parts generated at three loops by the unitarization prescription 
Eq. (4.9). The latter contributes to \M\ 2 at 0(a^ n e 4 ), i.e. it is suppressed in powers of 
a wn , but enhanced in e. By investigating the imaginary part of the dominant isospin 1 = 
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amplitude Aio(s), with Eq. (C.l) added appropriately to the representation Eq. (4.4), we 
find that the two-loop imaginary part is suppressed by more than a factor of 30 relative to 
the one-loop piece at the center of the Dalitz plot, and by roughly a factor of 2 relative to 
the three- loop part. This suppression grows even stronger when considering derivatives of 
the amplitude, as s-dependence in ImF(s) is even further suppressed in e, see Eq. (C.l). 
We therefore neglect these terms of 0(ia^ n e ) in our analysis, and consider their effects 
to be safely included in our error estimates due to partial higher-order resummation. The 
smallness of the imaginary parts due to three-pion cuts is in accordance with findings from 
ChPT at two loops [6] as well as from dispersion relations [22] . 
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